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Abstract. Wc prove the analytic existence of a symmetric periodic simulta- 
neous binary collision orbit in a regularized planar pairwisc symmetric equal 
mass four-body problem. Wc provide some analytic and numerical evidence 
for this periodic orbit to be linearly stable. We then use a continuation method 
to numerically find symmetric periodic simultaneous binary collision orbits in 
a regularized planar pairwise symmetric 1, m, 1, m four-body problem for m 
between and 1. We numerically investigate the linear stability of these peri- 
odic orbits through long-term integration of the regularized equations, showing 
that linear stability occurs when 0.538 < m < 1, and instability occurs when 
< m < 0.537 with spectral stability for m Ri 0.537. 



1. Introduction 

In the TV-body problem, linearly stable periodic orbits trap around themselves 
bounded, non-chaotic motion of the N masses [4]. Some of the known examples of 
linearly stable periodic orbits in the three-body problem are the elliptic Lagrangian 
triangular periodic orbits for certain values of eccentricity and the three masses 
[9], [15], and the Montgomery-Chenciner figure-eight periodic orbit for equal masses 
[3], [11], [16], [7], [8]. Other examples of linearly stable periodic orbits in the three 
or four-body problem involve binary collisions (BC) and/or simultaneous binary 
collisions (SBC). The regularization of these kinds of singularities is achieved by a 
generalized Levi-Civita type transformation and an appropriate scaling of time, as 
adapted from Aarseth and Zare [1] to the particular problem (also see [13]). 

Schubart [21] numerically discovered a singular symmetric periodic orbit in the 
collinear three-body equal mass problem. In this orbit, the inner mass alternates 
between binary collisions with the two outer masses. Henon [5] extended Schubart 's 
numerical investigations to the case of unequal masses. Only recently did Vcnturelli 
[24] and Moeckcl [10] prove the analytic existence of the Schubart orbit when the 
outer masses arc equal and the inner mass is arbitrary. The linear stability of the 
Schubart orbit was determined numerically by Hietarinta and Mikkola [6] revealing 
that linear stability occurs only for certain choices of the three masses. Numerically, 
non-Schubart-like linearly stable periodic orbits in the collinear three-body problem 
were found by Saito and Tanikawa for certain choices of the three masses [17], [18]. 

Sweatman [22], [23] and Sekiguchi and Tanikawa [19] numerically found a sym- 
metric Schubart-like orbit in the symmetric collinear four-body problem with masses 
1, m, m, and 1. This Schubart-like periodic orbit alternates between simultaneous 
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binary collisions of the two outer pairs of masses and binary collisions of the inner 
two masses. Ouyang and Yan [14] proved analytically the existence and symmetry 
of this orbit. In the regularized setting, this periodic orbit has a symmetry group 
isomorphic to D2 , of which both of the generators are time-reversing symmetries. 
(Here the dihedral group is the group of symmetries of the regular fc-gon.) 
Sweatman [23] numerically showed that the Schubart-like orbit is linearly stable 
when < m < 2.83 or m > 35.4, and is otherwise unstable. This linear stability 
was confirmed in [2] using Robert's symmetry reduction technique [16]. 

Ouyang, Yan, and Simmons [12] numerically found and analytically proved the 
existence and symmetry of a singular symmetric periodic orbit in the fully symmet- 
ric planar four-body problem with equal masses. (In the fully symmetric planar 
four-body equal mass problem, the position of one mass determines the positions 
of the other three masses.) In this orbit, the four masses alternate between differ- 
ent simultaneous binary collisions. In the regularized setting, this periodic orbit 
has a symmetry group isomorphic to D4, of which one of the two generators is a 
time-reversing symmetry, while the other generator is a time-preserving symmetry. 

In this paper we consider the existence and linear stability of time-reversible 
periodic simultaneous binary collision orbits in the planar pairwise symmetric four 
body problem. The positions of the four bodies in the plane are {xi,X2), {x3,x4), 
{~xi, —X2), and (— X3,— 2:4), where the corresponding masses are 1, m, 1, m with 
< TO < 1. With t as the time variable and ' = d/dt, the momenta for the four 
masses are (wi,W2) = '2{xi,X2), {103,014} = 2m{x3,X4), —(011,102), and — (w3,a;4). 
The Hamiltonian for the pairwise symmetric planar four-body problem is H = 
K — U, where 



and 



U= , „ = + 



4m ' 
2m 



2^yxl + xl \/{X3-Xiy + {X4-X2Y 

2m 



(Xi + X^Y + {X2 + X4Y 2sjx\ + X4 

The angular momentum for the pairwise symmetric planar four-body problem is 

A = X\L02 — X2OJ1 + X3U4 — X4UJ3. 

The center of mass is fixed at the origin, and the linear momentum is zero. With 

/ 
-/ 

for I the 4x4 identity matrix, the vector field for the pairwise symmetric pla- 
nar four-body problem is JVH, i.e., the Hamiltonian system of equations with 
Hamiltonian H are Xi = dH/duJi, uii = —dH/dxi, i = 1,2,3,4. 

The initial conditions for the orbits of interest has the first body of mass 1 located 
on the positive horizontal axis with its momcntiim perpendicular to the horizontal 
axis, and the first body of mass m located on the positive vertical axis with its 
momentum perpendicular to the vertical axis. Specifically, at f = we have 

x\ > 0, X2= 0, X3 = 0, X4> 0, with X4< Xl, 
coi =0, CO2 > 0, CO3 > 0, 0J4 = 0, with i02 < 0J3, 
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Figure 1. The symmetric periodic simultaneous binary collision 
orbit in the planar pairwise symmetric four-body problem for m = 
1 (left) and m = 0.539 (right). The two red curves are those traced 
out by ±{xi{t), X2{t)), and the two blue curves are those traced out 
hy±{x3it),Xiit)). 

at which H is defined. The first objective is to find, for < m < 1, values of 
xi, X4,ui2i'^3 a.t t — such that (i) x^ ~ xi = and X4 — X2 = with x^ -\- x\ ^ Q 
at some t = to > 0, (ii) Xi + x^ = and X2 + x^ — Q with x\ -\- x\ ^ Q ai some 
t = ti > to, (iii) the orbit extends to a symmetric periodic orbit, and (iv) the 
periodic orbit avoids all the other kinds of collisions. Such an orbit experiences a 
simultaneous binary collision in the first and third quadrant at t = to, and then 
another simultaneous binary collision in the second and fourth quadrants at t = ti, 
before returning to its initial conditions at some t — t2 > ti. The presence of 
collisions along the orbit necessarily imposes zero angular momentum on the orbit, 
thus requiring that X1UJ2 — X4LO3 = at t = 0. Examples of these symmetric periodic 
simultaneous binary collision orbits in the planar pairwise symmetric four body 1, 
TO, 1, TO problem are illustrated in Figure 1 for to = 1 and to = 0.539. The second 
objective is to investigate the linear stability of the symmetric periodic simultaneous 
binary collision orbits as to varies over interval (0, 1]. 

The regularization of the simultaneous binary collisions, as described by (i) and 
(ii) above, in the Hamiltonian system of equations with Hamiltonian H plays a 
key role in achieving the two objectives. Section 2 details this regularization which 
consists of two canonical transformations followed by a scaling of time t — 6{s) 
with s as the regularizing time variable, producing a new Hamiltonian T in the 
extended phase space. Section 3 describes a scaling of orbits of the Hamiltonian 
system of equations with Hamiltonian T which shows that any such periodic solution 
always belongs to a one-parameter family of periodic solutions for which the linear 
stability is the same for all periodic solutions in the family. Section 4 describes the 
symmetries of the Hamiltonian T which are used to construct periodic solutions 
with a D4 symmetry group generated by a time-reversing symmetry and a time- 
preserving symmetry. 

In Sections 5 and 6, we prove the analytic existence of a periodic simultaneous 
binary collision orbit 7(5), with a D4 symmetry group, for the Hamiltonian system 
of equations with Hamiltonian F with to = 1, and investigate its linear stability. 
The proof extends the analytic existence of a symmetric periodic simultaneous 
binary collision orbit in the fully symmetric planar four-body equal mass problem, 
as found in [12], to the planar pairwise symmetric four-body equal mass problem. 
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The symmetric periodic simultaneous binary collision orbit in the fully symmetric 
four-body equal mass problem is known to be linearly stable [2], which provides 
some analytic evidence for the linear stability of 7(s). We give numerical evidence 
that supports the linear stability of 7(5). 

In Section 7, we numerically continue the symmetric periodic simultaneous bi- 
nary collision orbit 7(5) from to = 1 to a symmetric periodic simultaneous binary 
collision orbit j{s; to) for to < 1, and then investigate the linear stability of 7(5; to) 
as TO varies in the interval (0,1]. Wc use trigonometric polynomials as approxi- 
mations of the periodic orbits (cf . [20] ) . The numerical algorithm for continuation 
starts with a trigonometric polynomial approximation of j{s; 1) = 7(5) that is used 
as a initial guess for a trigonometric polynomial approximation of 7(5; 0.99), where 
the coefficients of the trigonometric polynomial are optimized through a variational 
approach. This process is repeated, using the optimized approximation of 7(5; 0.99) 
as the initial guess for 7(5; 0.98), etc., until an optimized approximation of 7(5; 0.01) 
is obtained. Numerical integrations of the Hamiltonian system of equations with 
Hamiltonian T show that 7(5; to) is linearly stable when 0.54 < m < 1, and is 
unstable when < to < 0.53. Wc then conduct a refined continuation for values 
of TO between 0.53 and 0.54 at 0.001 increments. Numerical integration of j{s; to) 
for these values of m reveals that linear stability appears to hold for 7(5; to) when 
0.538 < TO < 1 and that instability appears to hold when < to < 0.537 with 
spectral stability for to « 0.537. 

2. Regularization 

We adapt the regularization of Aarseth and Zarc [1] to the planar pairwise sym- 
metric four-body problem to regularize simultaneous binary collisions as described 
in the first objective. The first canonical transformation is 

{xi,X2,X3,X4,uJi,0J2,uJ3,uJ4) {qi, 92, 93, 94, hi, h2, hs, Ha) 
determined by the generating function 

Fi{xi, X2, X3, Xi, hi,h2, ha, /I4) = hi{xi-X3) + h2{X2-X4) + h3{xi+X3) + h4{X2+X4). 

So the first canonical transformation is 

wi = - — =hi+ hs, UJ2 = ^ — = h2 + hi, 

ox\ 0x2 

= -7, — = -hi + hs, Ui = - — = -h2 + hi, 

axs oxi 

dFi dFi 

91 = =Xi- X3, 92 = -WT- =X2- Xi, 

ohi oh2 

dFi dFi 

93 = = Xl + X3, 9i = = X2 + Xi. 

ohz ohi 
91 +93 92+ 94 93 -91 94- 92 

Xl = , X2= , X3 = , Xi= , 

Wl - W3 UJ2- LOi UJi+U)3 W2 + W4 

hi = , n,2 = , 113 = , hi = . 

^2 , ^2 _ (gi + 93)^ + {92 + 94)^ ^2,9_ {93 -91? + {94 -92? 

Xi-\- X2 — , X3-\-Xi — - . 



Here 

and 

So 
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The new Hamiltonian is H = K — U, where 



^ {hi + h^f + ih2 + h^f ^ (/l3 - hif + (/14 - h2f 

K = 1 



4m 



and 

U = 



2m 



2m 



Vi9l + 93? + (52 + 54)2 ^g2 ^ g2 ^g2 ^ ^2 ^(gi - 33)2 + - 54)2 " 

The second canonical transformation is 

{9l,92,g3,94,hi,h2,h3,h4) (ui,U2,M3,'"4,'?^1,^^2, V3,i;4) 

determined by the generating function 

4 

F2{hi,h2,h3,h4,Ui,U2,U3,U4) = y^^hjfj{ui,U2,U3,U4), 

i=i 

where 

/l = _ uj, f2 = 2uiW2, /a = "3 - ^4) /4 = 2U3W4. 

So the second canonical transformation is 
dF2 
dhi 

dF2 
dui 
dF2 



91 

93 
Vi 

V3 



2 2 
Ui - U2, 



U3 - W4> 



2hiUi + 2/12M2, 
2/13^3 + 2/14U4, 





dF2 


52 = 


dh2 




dF2 


54 = 


dh4 




dF2 


V2 = 


dU2 




dF2 


V4 = 


dU4 



2uiU2, 

2U3U4, 

—2h\U2 + 2/12W1, 
—2hsU4 + 2/14^3. 



Here 
and 



5? + 52 = (w? - ujf + A.u\ul = {u\ + ulf, 



2,2 I 2 2\2 I A 2 2 / 2 I 2\2 
93+94 = («3 - "4) + 4^3% = + "4) • 

Also, solving §2 = 2uiU2 for U2 and substituting this into gi=u\— gives 

«f-5i-?-f =0. 

Solving this quadratic in u\ and noting that > gives 

2 _ 51 + ^51 +52 
^ • 

Substituting this into gi=u\ — and solving for gives 

, 2 _ -51 + V5i + 52 



Similarly 



53 + \/93 + 54 
2 



and 



-53 + v'5i + 54 
2 



The equations in fi, W2, ^3, W4 arc linear in /ii, /i2j ^3, ^4: 





= 2 


Ml 


M2 




'hi 




V3 


= 2 


W3 


1*4 




'hs 






-^2 


Ml 




>2. 








— U4 


U3 




/14 
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Solving these equations for fti, /i2, /14 gives 





'hi 




1 


Ui 


-U2 




Vl 


1 


viUi — V2U2 






2{ul + ul) 


U2 


Ui 






~ 2{ul+ul) 


V1U2 + V2U1 


and 






1 


U3 


— U4 






1 


V3U3 - V4U4 




hi 


2(ui + «2) 




U3 










The new 


Hamiltonian is H - 


= k- 


U, where 






16 V 




'{vf +v 






) 


+ ( 


vi + vi){ui + uiy 




m J 















1 / (^3^3 - V4U4){viUi - V2U2) + <yV3U4 + V4U3){viU2 + V2U1) 

SV m) {ul + ul){ul + ul) 



and 



U = 



1 



+ 



2m 



+ 



2m 



-ui+ui- uiY + (2uiM2 + 2U3U4Y ul + ul ul + ul 

2 

m 

+ 



{ul - u| - U§ + + (2uiU2 - 2U3M4)2 

Introduce a new time variable s by the regularizing change of time 



dt 
ds 



{ui + ul){ul + ul). 



To simplify notation, set 



Ml = viui — V2U2, 
M3 = V3U3 - V4U4, 

= Ul — U2 + u'^ — U4, 

Mr = Ml — u| — U3 + U4, 



M2 = V1U2 + V2U1, 
M4 = V3U4 + V4U3, 
Me = 2uiU2 + 2U3U4, 
Mg = 2uiU2 — 2U3U4. 



The Hamiltonian in the extended phase space with coordinates Ui, U2, U3, U4, E, 
Vl, V2, V3, V4, t is 

f = £ - ^) =1 (1 + 1) + vDiul + ul) + {vl + vl){ul + ul) 

I (M3M1 + M4M2) 



+ 



1 - 



{ul+ul){ul + ul) ,2 +^2,2 2. 

^ M'i + M| 2mK+U2 + «3+%) 
m'^{ul+ul){ul+ui) 
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With ' = d/ds, the Hamiltonian system of equations with Hamiltonian f is 

+ (M3U1+M4U2), 

- ) V2{ul + ul) + \(l- - ) (M3(-U2) + M4U1), 





_ dt 

dvi 


1 

~8 


\ m 




dt 


1 

~8 


(i + i 

\ m 


U3 


_ dt 

dv3 


1 

^8 


+ - 

\ m 


U4 


_ dt 

dV4 


1 

~8 


\ m 


and 








v[ = 


dt 

dui 







vsiul + ul) + ^(^l- ^ j (Ml Us + M2W4), 

V4{ul + ul) + l (1- — ) (Mi(-U4) + M2M3), 

by ui J 



+ 



2ui(m| + uf) _ 2('(if + + ul){M5Ui + MQU2) 

^ ^/M'i + Mi [M| + M|]3/2 

2m2ui(ui + ul) _ 2m2(uf + u2)(u2 + ^^(Mytti + M8U2) 
^M^ + Mi [M2 + M|]3/2 

+ 4mMi + 2£;ui(u3 + W4), 

= = - ^ (1 + 1) U2(.3^ + .1) - ^ (1 - ^) {Ms{-V2) + M4.1) 

2m2(M3+M4) _ 2(Mf + ui)(M§ + uI){M5{-U2) + Mgm) 

VM| + M| [M| + M|]3/2 

2m'^U2{ul + ul) _ 2Tr?{u\ + + «1)(M7(-U2) + MgUi) 

+ 4mu2 + 2Eu2{ul + uj), 
, __dt 1 1 \ . , o. 1 1 



'3 



1 + + ^^i) - ^ (^1 - ^) (^1^^3 + M2V4) 

^ 2u3{ul + ul) 2{ul + + ul) + MeUi) 



+ 



VMi + Mi [M| + M|]3/2 

+ M| [M2 + M2]3/2 

+ 4mu3 + 2^W3(u? + W2), 

2^4 (m^ _ 2{ul + ul){ul + ul){M5{-U4) + MeUj) 

^ VMi + Mi [M| + M|]3/2 

2m'^Ui{ul + ul) 2m?{ul + ul){ul + ul){MrU4 + Ms{-U3)) 
^ ^yM^+Ml [M| + M|]3/2 

+ 4mu4 + 2Eu4iul +ul), 
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along with the auxihary equations, 

E'=|^=0, t' = -^ = {u\+ul){ul+ul). 

On the level set T = 0, the value of the Haniiltonian H (i.e., the energy) along 
solutions of the Hamiltonian system of equations with Hamiltonian F is E. Indepen- 
dent of the values of E and F, the angular momentum A = xiOJ2—X2i^i+X3U)4—XiUi3 
in the coordinates ui, M2, U3, U4, wi, t)2, ^3, V4 simplifies to 

vl = i [ - V1U2 + V2U1 - V3U4, + V4,U3] . 

On the level set F = 0, two simultaneous binary collisions in the planar pairwise 

symmetric four-body problem have been rcgiilarizcd in the Hamiltonian system of 
equations with Hamiltonian F. The simultaneous binary collision — xi =0 and 
X4 — X2 = with x1+ X2^ corresponds to wf + tt| = with u'^ + 0. These 
imply that M| + M| = 4(a;f + x^) and + = A{xl + xl) ^ 0. From 
f = it follows that 

Similarly, the simultaneous binary collision 2:3 + a; 1 =0 and X4 + X2 = with 
xf + xl ^ corresponds to u^+ul — with + 7^ 0, and hence that M^ + Mq = 
4{xl + xl) ^ 0, M| + M| = 4(x? + xl) ^ 0, and, from f = 0, that 

2 2 32m2 

V3 +vi- 



m+1 

On the level set F = 0, the other singularities of the planar pairwise symmetric 
four-body problem have not been regularized in the Hamiltonian system of equa- 
tions with Hamiltonian F. The binary collision Xi = 0, X2 = with .T3 + a;| ^ 
corresponds to M| + = and M| -h M| ^ with + m| ^ and uj + ul^ 0. 
The binary collision 2:3 = 0, a;4 = with -I- x| 7^ corresponds to M| + Mg 7^ 
and M7 + M| = with + ul and ul + ul ^ 0. Because of the pairwise 
symmetry, there are no triple collisions. Total collapse xi — 0, 2:2 = 0, 2:3 = 0, 
X4 = corresponds to ui = 0, U2 = 0, 1x3 = 0, and U4 = 0. A solution of the 
Hamiltonian system of equations with Hamiltonian F is called nonsingular if it 
avoids the unregularized binary collisions and total collapse singularities. 

We establish next the correspondence between the original coordinates and the 
regularized coordinates for the initial conditions given in the Introduction. 

Lemma 2.1. The conditions {att = 0) 

Xl > 0, a;2 = 0, 2:3 = 0, 2:4 > 0, with 2:4 < 2:1, 
LOi =0, UJ2 > 0, C03 > 0, W4 = 0, with 0J2 < 

correspond to the conditions {at s = 0) 

U3 = ±ui, U4 = with U1U2 < 0, \u2\ < (\/2 — l)|ui|, 
vs = TVi, V4 = ±V2, with < V1U2 + V2U1 < V2U2 — viui. 
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Proof. Suppose that X2 = 0, X3 = 0, wi = 0, ^04 = 0, and that xi, X4, u)2, and 
are positive with X4 < Xi and W2 < ^3. The first canonical transformation imphes 
that 

9i= 93> 0, 92 = -94 < 0, with I52I < ffi, 
hi = < 0, h2 = h4 > 0, with ft2 < l^-il- 

The second canonical transformation can be rendered in the following complex 
notation identities: 

gi + ig2 =u\-u% + 2iuiU2 = {ui + iu2Y, 
93 - i9i =u\-u\- 2iu3U4 = (us - 1^4)^, 

_V2U2 — ViUi .V1U2+V2U1 

{—Vi + iv2){ui — iu2) —Vi + iv2 



and 



2{ul + ul) 2(ui+m2)' 



_ V3U3 - V4U4 . V3U4 + V4U3 
2{ul + ul) 2(^2 + ^2) 

{V3 + iV4){u3 + iU4) V3+iV4 



2(u§ + ul) 2(u3 - iu4) ' 

The two identities relating 91,92,93,94 with ^1,^2,^*3,^4 imply that 

{ui + iu2f =91+ i92 =93- i9i = {u3 - iu4f. 
Thus Ul + iu2 = ±(w3 — iu4), and so 

Us = ±Ui, U4 = TU2- 

Now 52 = 2uiU2 and 52 < imply that 

U1U2 < 0. 

For a complex number z, let arg(z) denote the argument of z, i.e., the angle z 
makes with the positive horizontal axis, modulo 2tt. Since gi > 0, 92 < 0, and 
1521 < 9i, then arg(5'i + 1172) € [77r/4, 27r). Since {ui + iu2)^ = <?i + 192, it follows 
that arg(wi + iu2) e [Ttt/S, tt) U [IStt/S, 27r). This implies that 

1^2! < tan(7r/8)|ui|. 

It is easily shown that tan(7r/8) = \/2 — 1. 

The equalities hi = — /13 and h2 = /14 imply that —hi+ i/12 = /13 + i/14. By the 
identities relating /ii, /i2, /is, /i4 and ui,U2,U3,U4,vi,V2,V3,V4, we have that 

-Wl + iV2 _ V3+ iV4 

2{ui + iu2) 2{u3 - iu4) ' 
Since wi + iu2 = ±(u3 — iu4), then —vi + iv2 = ±(f3 + W4), and so 

Since h2 < \hi\, hi < 0, and ft-2 > 0, the second canonical transformation implies 

f 1M2 + f 2tti _ _ f 1^2 - t^lM2 

^ 2(^x2 + ^2) -''^S fti- 2(n2+u2)- 
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From this it follows 

< V1U2 + V2U1 < V2U2 — viui. 

Now suppose that U3 — ±ui, U4 = TU2, U1U2 < 0, \u2\ < (\/2 — 1)I"i|j ^3 — T^i, 
V4 = ±V2, and < V1U2+V2U1 < V2U2—V1U1. The second canonical transformation 
implies that 

gi = ul- ul, g2 = 2uiU2 < 0, 53 = - ul, 54 = -2uiU2 > 0. 

The first canonical transformation implies 

xi = — U2 > 0, X2 =0, X3 = 0, X4, = —2uiU2 > 0. 

In complex notation, xi — ix4 = (ui +^2)^- Then \u2\ < (\/2— l)|'Ui| implies that 
\x4\ < tan(7r/4)|a;i| = |a;i|, i.e., 

X4 < Xi. 

The equalities Us = ±ui, U4 = =FW2, = T^^i, and V4 = ±V2 imply that uf + ^2 = 
W3 + ul, and 

V3U4 + V4US = V1U2 + V2U1, V3W3 — U4M4 = —ViUi + V2U2- 

Thus the second canonical transformation implies 

_ _ ~VlUi + ^2^2 _ V3.U3 - V4U4 _ 

_ V1U2 + V2U1 _ V3U4 + V4U3 _ 
' ~ 2{ul + ul) ~ 2{ul + ul) ~ 
Furthermore, the inequalities < V1U2 + V2U1 < V2U2 — f iWi imply that 

-hi >h2>0. 
The first canonical transformation now implies that 

cji = hi + hs = 0, CJ2 = h2 + h4> 0, W3 = —hi + /13 > 0, W4 = — /12 + /14 = 
with W3 = —2hi > 2/12 = 002- n 

3. A Scaling of Periodic Orbits and Linear Stability 

A certain scaling of solutions of the Hamiltonian system of equations with Hamil- 
tonian T produces more solutions. When applied to a periodic solution, this scaling 
leads to a one-parameter family of periodic solutions. The proof of the following 
result is a straight-forward verification. 

Lemma 3.1. If j{s) = (ui(s), 1*2(5), M3(s), ^4(5), 1^1(5), ^2(5), t'3(s), i'4(s)) is ape- 
riodic solution of the Hamiltonian system of equations with Hamiltonian T on the 
level set F = with period T and energy E, then for every e > 0, the function 

7e(s) = (eui(es), eM2(es), eu3(es), eu4{es), wi(es), W2(es), V3{es),V4{es)) 

is a periodic solution of the Hamiltonian system of equations with Hamiltonian F 
on the level setT = Q with period = e~^T and energy = e~'^E. 

The linear stability of a periodic orbit of the Hamiltonian system of equations 
with Hamiltonian F is determined by the linearization of the equations along the 
periodic orbit. By Lemma 3.1, a periodic orbit 7(3) on the level set F = with 
period T and energy E embeds into a one-parameter family 7e(s) of periodic orbits 
on the level set F = with period and energy E^. The linearization of the 
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Hamiltonian system of equations with Hamiltonian f along the periodic orbit 7e(s) 
is 

X' = JV2f(7e(s))X 

where V^F is the matrix of second-order partials of F. Let Xg{s) be the solution 
of the linearization of the equations along je{s) that satisfies ^e(O) = / (the 8x8 
identity matrix). The nionodromy matrix for 7e(s) is Xg{T^), and the eigenvalues 
of Xg{T^) are the characteristic multipliers of 7e(s). A characteristic multiplier 
A of 7e(s) is defective if its geometric multiplicity is smaller than its algebraic 
multiplicity, i.e., its generalized eigenspace Uj>iker(Xe(Te) — Xiy is not the same 
as its eigenspace keT{Xf^{T^) — XI). 

Lemma 3.2. If"f{s) is aperiodic orbit of the Hamiltonian system of equations with 
Hamiltonian T on the level set F = 0, then for each e > 0, the periodic orbit 7c (s) 
has 1 as a defective characteristic multiplier with algebraic multiplicity at least two. 

Proof. By Lemma 3.1, 7e(s) is a one-parameter family of periodic orbits on the level 

set r = with period = e~^T and energy = ^^^E, where 71 (s) = 7(5) and 
T is the period and E is the energy of 7(3). For each e > 0, the periodic solution 
7e(s) satisfies 7j'(s) = JVT{'y^{s)), and so 

and 

^7.) (s) = JV2f(7.(.))(^|7e)(^). 

Since X^ is a fundamental matrix solution that satisfies ^e(O) — I, it follows that 

d 



jUs) = X,{sWM and ^7.(5) = Ms) ( 



) (0)- 



The T^-periodicity of 7e(s) implies the T^-periodicity of 7^(s), and so 

7^(0)=7^(T,) = X,(T,)7^(0) 

where 7^(0) 7^ 0. Thus 1 is an characteristic multiplier of 7e(s). Since Xcf Tf ) is 
symplectic, the algebraic multiplicity of this characteristic multiplier is at least two. 
Now 

7j'(s) = (e^Ui(es), . . . , e^U4(es), eu'(es), . . . ,eW4(es)) 

and 

d 



Q^lej (s) = {ui{es) + esu[{es),...,U4,{es)+esu'4{es),sv[{es),...,sv'^{es)). 
The T^-periodicity of 7e(s) and 7j'(s) implies that 

Thus 

Xe(T.) (|7e) (0) = (|7.) (T.) = (|7.) (0) + e-T.7;(0). 
Consequently, 

(X,(T,) - I) (^7e) (0) = e-^T,7,'(0) ^ 
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and 



(X,(T,)-7r (^7.)(0) = 0. 



This shows that 

U ker(X,(T,) - ly ^ ker(X,(T,) - 7), 
i.e., that 1 is a defective characteristic multipUer of 7e(s) for each e > 0. 



□ 



For each e > 0, the periodic orbit 7e(s) is spectrally stable if all of its character- 
istic multipliers have modulus one. By Lemma 3.2, the periodic orbit 7e(s) has 1 as 
a defective characteristic multiplier with algebraic multiplicity at least two, and so 
the monodromy matrix Xf{Tg) it not semisimple. However, as shown in the proof 
of Lemma 3.2, the two-dimensional subspace 

Span (jM, ^ ^ 



de 



7. (0) 



is X(;(Tg)-invariant. The periodic orbit 7e(s) is said to be linearly stable if it is 
spectrally stable and there exists a 6-dimensional Xe(Te)-invariant subspace U2 
such that Ui + U2 = and X^(T^) restricted to U2 is semisimple. 

Theorem 3.3. Suppose 7(5) is a periodic orbit of the Hamiltonian system of equa- 
tions with Hamiltonian T on the level set F = 0. Then 7e(s) is spectrally {linearly) 
stable for some e> if and only i/7e(s) is spectrally {linearly) stable for all e > 0. 

Proof. Let 7(5) be a periodic orbit with period T and energy E. For each e > 0, 
the periodic orbit 7e(s) has period = e~^T and energy = e~'^E. Note that 
71 (s) = 7(s), Ti = T, and Ei = E. 

It suffices to show that X^{Tf) is similar to Xi{Ti) for all e > 0. The fundamental 
matrix X^{s) satisfies X'^{s) = \/'^r{je{s))X^{s), X^{0) = I. Since the matrix V^f 
is symmetric, there are 4x4 matrix functions A,B,D that satisfy 

A Si 



D 



with A"^ 
A = 



= A and 

duidu-i 



D. Here 

, B = 



i,j=l,...,4 



duidvi 



D = 



dvidv-i 



J <,j=l,...,4 



From these, it is straightforward to show that 

A(7,(s)) = A(7i(es)), B(7e(s)) = eB(7iM), D{^,{s)) = e^D{^,{es)), 

where the value of energy on the left of each of these is E^ and on the right is Ei . 
Thus 

'eB^{ji{es)) e^D{j,{es))- 
_-A(7i(es)) -65(71 (es))_ " 
For I the 4x4 identity matrix, define the nonsingular matrix 

-1/2/ 



JV2f(7.(s)) = 







Then 



eF-Wf(7i(es))Fe = e 



■ i?^(7iM) eD{^,{es)y 
-e-iA(7i(es)) -B(7i(es)). 



JV2f(7e(5)). 
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Thus 

{Yf'Xi{e.s)Y,y = er-i JV2f(7i(e.s))Xi(es)r, 

= eY-^JV^r{ji{es))Y,Y-'X,{es)Y, 

= JV2f(7,(s))y-iXi(es)y,. 

Since Xt{s) and Y~^Xi{es)Y^ both evaluate to the 8x8 identity matrix at s = 0, 
uniqueness of solutions implies that Xf{s) = Y~^Xi{es)Yf for all s. In particular, 

X,{T,) = Yf^Xi{eT,)Y, = Yf'X,{T^)Y,. 

Therefore, Xe(Te) and Xi{Ti) are similar. □ 

4. Symmetries 

The Hamiltonian system of equations with Hamiltonian F has a group of sym- 
metries isomorphic to the dihedral group D4 

-1 



64 = (a6)^ = e). With 







G = 



define the matrices 






F 








-F 




















F 








-F 






Sr 



1 
1 

-G 








G 







G 







-G 



These matrices satisfy Sp = -I, Sp = I, Sq = I, and (SfSg)'^ = I- Fixing the 

value of E, these matrices satisfy F o S'p = F and F o 5^ = f, and so Sp and Sq 
are the generators of the D4-symmetry group for F. If 

lis) = {ui{s),U2{s),U:i{s),U4{s),Vi{s),V2{s),V3{s),V4(s)) 

is a solution of the Hamiltonian system of equations with Hamiltonian F, then 
5f7(s), 51,7(5), and Sgj{—s) are also solutions of the Hamiltonian system of 
equations with Hamiltonian F. This means that Sp is a time-preserving symmetry 

and that Sq is a time-reversing symmetry. 

Lemma 4.1. If for some sq > there is a nonsingular solution j{s), s S [0, sq], 
of the Hamiltonian system of equations with Hamiltonian T such that for constants 
Ci 7^ 0; C2 7^ 0, pi ^ 0, and /O2 7^ there holds 

ui{0) = Ci, U2{0) = C2, U3(0) = Ci, MO) = -C2, 

^^l(0) = Pi, V2{0) = p2, Vs{0) = -pi, V4{0) = P2, 

and 

ui(so) = 0, U2(so) = 0, U3(so) ^ 0, M4(so) 7^ 0, 

vi{so) 0, W2(so) 7^ 0, V3{so) = 0, 1^4(50) = 0, 

then 7(5) extends to a periodic orbit with period 8sq and a symmetry group isomor- 
phic to D4 such that 

wi(3so) 7^ 0, U2(3so) 7^ 0, U3(3so) = 0, U4(3so) = 0, 
vi{3so) = 0, V2{3so) = 0, t;3(3so) 7^ 0, Vi{3so) 0, 
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and 

ui(5so) = 0, U2(5.so) = 0, ?i3(5so) ^ 0, U4(5so) ^ 0, 
vi{5so) ^ 0, v-2.{^sq) ^ 0, W3(5so) = 0, -^4(550) = 0, 

and 

Wi(7so) 7^ 0, -^2(750) 7^ 0, ?i3(7so) = 0, ^4(730) = 0, 
t;i(7so) = 0, U2(7so) = 0, t;3(7so) ^ 0, t;4(7so) ^ 0. 

Proof. The curve 5g7(2so — s), s € [sq, 2so], is a nonsingular solution of the Hamil- 
tonian system of equations with Hamiltonian F whose conditions at s = Sq are 

5,37(250 - So) = Sgi{sq) = (0,0,U3(so),U4(so),ui(so),t;2(so),0,0) = 7(so). 

By uniqueness of solutions, 5g7(2so — sq), s e [so,2so], is the extension of 7(s), 
s e [0,so], to [so,2so], i.e., for s e [so,2so], 

ui{s) = — ui(2so - s), U2{s) = -M2(2so - s), 

Us{s) = U3(2so - S), U4{s) = 1t4(2.So - s), 

vi{s) = Vi{2so - s), V2{s) = V2{2so - s), 

^3(5) = -W3(2so - s), Vi{s) = -V4{2so - s). 

The conditions of this extension at s = 2so are 

5g7(2so -2so) = 507(0) = (-Ci, -C2, Ci, -C2,Pi,P2, Pi, -^2)- 

The curve Sp'yis — 2so), s G [2so,4so], is a nonsingular solution of the Hamil- 
tonian system of equations with Hamiltonian F. Its conditions at s = 2so are 

S'f7(2so - 2so) = S'f7(0) = (-Ci, -C2, Ci, -C2,Pi, P2,Pi, -^2)- 

By uniqueness of solutions, S'f7(s — 2so), s £ [2so,4so], is the extension of 7(5), 
s e [0,2so], to [2so,4so], i.e., for s G [2so,4so], 

ui{s) = -U3{s - 2sq), U2{s) = U4{s - 2so), 
U3{s) = ui{s - 2so), Ui{s) = -U2{s - 2so), 

vi{s) = -v-i{s - 2so), V2{s) = V4{s - 2so), 
vsis) = vi{s - 2so), V4{s) = -V2{s - 2so). 

The conditions of this extension at s = 4so are 

5ir7(4so - 2so) = Sfj{2so) = SpSc^iO) 

= (-Cl, -C2, -Cl> C2, -Pl, -P2,Pl, -P2). 

The curve Sp'y{s — 4so), s G [4so,8so], is a nonsingular solution of the Hamil- 
tonian system of equations with Hamiltonian F. Its conditions at s = 4so are 

51,(430 - 4so) = 5|7(0) = (-Ci,-C2,-Ci,C2,-pi,-P2,Pi,-p2). 



PERIODIC SBC ORBITS IN THE PLANAR PAIRWISE SYMMETRIC PROBLEM 



15 



By uniqueness of solutions, Sp'y{s — 4so), s G [4so,8so], is an extension of 7(s), 
s e [0,4so], to [4so,8so], i-e., for s e [4so,8so], 

Ul{s) = -Ui{s - Asq), U2{s) = -U2{S - Asq), 

uais) = -U3{s - 4so), U4is) = -U4{s - 4so), 
wi(s) = -vi{s - 4so), V2{s) = -V2{s - 4so), 
n{s) = -^3(5 - 4so), ^4(5) = -V4{s - 4so). 
The conditions of this extension at s — Ssq arc 

SIj{8so - 4so) = SIj{4so) = S'|S'f7(2so) 
= {(1,(2X1, -(2, Pi, P2, -Pi, P2) =7(0)- 

Thus the overall extension of 7(s) from [0, sq] to [0, 8so] is nonsingular and periodic 
with period Ssq. By the construction of this extension, 

7(2so — s) = Sc'yis) and 7(5 + 2so) = Sp'jis) for all s, 

and so the periodic extension of 7(3) has a symmetry group isomorphic to D4. 
Simultaneous binary collision occur at s = So by hypothesis, at s = 3so where 

7(3so) = Sf7{so) = (-W3(0), «4(0), 0, 0, 0, 0, ui(0), -i;2(0)), 

at s = 5so where 

7(5.so) = 52 7(so) = (0, 0, -U3(0), -^^4(0), -^;i(0), -^;2(0), 0, 0), 

and at s = Tsq where 

7(7so) = S|7(3so) = (u3(0),-U4(0),0,0,0,0,-i'i(0),^;2(0)). 

Therefore, the symmetric periodic extension of 7(5) has four distinct simultaneous 
binary collisions between s = and s = 8so- D 

5. Analytic Existence in the Equal Mass Case 

When m = 1, there is an additional symmetry in the positions of the four masses 
that reduces the planar pairwise symmetric four-body equal mass problem to the 
fully symmetric planar four-body equal mass problem. We exploit this reduction 
to prove the existence of a symmetric periodic simultaneous binary collision orbit 
in the equal mass case. 

The additional symmetry is the Ansatz, 

X4 = Xi, X3 = X2, with \X2\ < Xi. 

From this it follows that 

UJ4 = uji, 0J3 = LU2 xi — X2 > 0, Xi + X2 > 0. 
From the first canonical transformation, we have 

gi=Xi- X2, g2=X2- Xi, 

g3=Xl+X2, g4^Xi+X2, 

, — L02 , a;2 — wi 

h, = h2 = , 

, Wi + W2 , Wi + W2 

h3 = , h4 = . 
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Prom the second canonical transformation, we have 



Thus 



2 1 + \/2 , . + -s/S / N 

Ul = ^ {X1-X2), U2 = [X1~X2), 

2 1 + a/2/ , \ 2 -1 + a/2/ , \ 

U3 = {X1+X2), M4 = {X1+X2). 



2uj 2u\ 2m? 



3 - "4 

x\ — X2 = 1= = a;i + a;2 — 



l + \/2 -l + \/2' l + \/2 -l + \/2" 

Since 2u\U2 = 92 = X2 — x\ <^ and 2U3M4 = 54 = Xi + a;2 > 0, it follows that 



U4 



U3 = (V2 - Ijws. 



l + \/2 

From the second canonical transformation, we have 

V\ = \/2(wi - U)2)ui, V2 = -(2 - V2){u)i - CJ2)ui, 

V3 = V2{u>i + UJ2)U3, W4 = (2 - V^)(wi + U)2)us. 



These imply that 

^ = (a;r - a;2)ni = ^^^^y ^ = (-^ + -^)-3 = 

and thus 

1^2 = ^7^^^ ^ -{V2 - l)vi, V4 = ^ V3 = {V2 - 1)V3. 

Substitution into the Hamiltonian system of equations with Hamiltonian F (and 
with m = 1) gives 

, 4-2^2 2 , (x/2-l)(4-2x^) 2 

Ml = ^ U2= t^lMs, 

; 4-2^ 2 / (y2-l)(4-2y2) 2 

% = ^ ^3%, W4 = ^ ^3%, 
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E' = 0, and 

(4 - 2y/2)uivl 



+ 4ui + 



4uiw| 



\A4+^ {u\ + utf^ 



4 



(72-1) 



(4 - 2V2)uivl ^ ^ 4ui«i 



+ 2(4 - 2V2)Euiul, 



4ufu§ 



3/2 



+ 2(4 - 2y/2)Euiul 

V "l + ^3 (WJ + U3 j 



(V2-1) 



(4 - 2^/2)u3^;? 



+ 4U3 



4UIU3 



4ufu| 



+ 2{A-2y/2)Eulu:i , 



t' ={4- 2^f2fu\u\. 

Because U2 = — (\/2 — l)wi, ?i4 = (v^— 1)m3, W2 = — (\/2 — and W4 = (v^— 1)113, 
the equations in v/'j, 1/4 , (,'2: ™d duplicate those in u'^, M3, v^, and U3. The Ansatz 
2:4 = xi, 0:3 = ./■2 with 1x21 < xi, therefore leads to the reduced system of equations, 



Ml 



4-2^2 



-■?^iW3> 



u, =- 



E' =0, 



4- 2v^ 



W3W1, 



(4-272)mi;| ^ ^ 4Mmi 
4 ' v^^^i^ 



4uft 



+ 2(4 - 2-^)Euyu\, 



, _ (4 - 2s/2)mvl 

Vo — — 

i 4 

i! =(4 - 2-^fulul. 

Scale the value of E by 



+ 4U3 + 



4MfU3 



3-^+2(4-2y2)^«>3, 



E 



E 



4- 2\/2' 



and define 



4- 2\/2, 



2u\u\ 



(4 - 2\/2)2£;u?u^. 



It is straight-forward to check that the reduced system of equations satisfies 



dr 

dvi 



Vi = - 



dui 



i = l,2, 
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and 

dt' dE 
Thus the system of reduced equations is HamiUonian. 

We win simpUfy the Hamiltonian F by a hnear symplectic transformation with 
a multipUer ^ 1. Define new coordinates {Qi, Q2, E, Pi, P2,t) by 

Qi Pi 



Q2 P2 

^3 = 7?ni^ ^3 = 



E= t = 2^/4(^/2- 1)3/2 r. 

(4 - 2^2)2 ^ ' 

This is a linear symplectic change of coordinates with multiplier 

1 

^~ 25/V\/2- 1' 

Under this linear symplectic transformation and the accompanying scaling a = s/ fj, 
of the independent variable s, the Hamiltonian T becomes 

r = ^ {P?Ql + PiQl) - ^/2(Q? + Ql) - jggl^ - EQlQl 

The reduced system of equations is the Hamiltonian system of equations with 
Hamiltonian F, 



da 



\P2Ql 



— ^ A Qr + 2^/2Q, + - ^^^-^ + 2i;g,Q„ 

— - -- A g. + 2^/2Q. + -^== - j^^^ + 2EQ2Q,, 

along with the auxiliary equations, 

The function F is a regularized Hamiltonian for the fully symmetric planar four- 
body equal mass problem with the bodies located at {xi,X2), (x2,xi), (— xi, —X2), 
and {—X2,—xi) (sec [12]). On the level set F = 0, the solutions have energy E. 
One regularized simultaneous binary collision occurs when Qi = and Q2 0, 
for which F = implies Pf = 16\/2, and for which the transformation between 
Qi,Q2 and xi,X2 implies xi — X2 = and xi + X2 7^ 0. The other regularized 
simultaneous binary collision occurs when ^ and Q2 = 0, for which F = 
implies P| = 16\/2, and for which the transformation between Qi,Q2 and xi,X2 
implies xi — X2 7^ and xi + X2 — 0. Total collapse occurs when Qi = and 
Q2 = 0, and is the only singularity in F that is not regularized. A solution Qi{a), 
Q2{cr), Piio'), ^2(0'), cr € [0,cro], for ao > 0, of the Hamiltonian system equations 
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with Hamiltonian T is nonsingular if it avoids total collapse, i.e., Qf + Q2 7^ for 
all a € [0,(To]. 

The following result is from [12]. The proof of it is a consequence of four equal 
mass bodies starting at (xi,X2), {x2,xi), {—Xi,—X2), and {—X2,—xi) with xi = 1 
and X2 = 0, and with the momenta (0, t?), (t?, 0), (0, — t?), and (—■(?, 0) for any d > 0, 
always having a simultaneous binary collision on the line X2 = xi at a time to > 
continuously depending on d, such that the cluster velocity ii(to) + ^2(^0) is a 
continuous function of 

Lemma 5.1. There exists i9 > 0, ctq > 0, and a nonsingular solution Qi{(T), Q2{cr), 
Pi{(j), P2{(^), cr G [0,(To]; of the Hamiltonian system of equations with Hamiltonian 
r on the level set F = such that 

Qi(o) = i, g2(o) = i, Pi(o) = -z?, P2(o) = i?, 

E=^-^^^^<0, r{o) = f^Q\{y)Ql{y)dy, 

and 

Qi(ao)=0, g2M>0, Pi (ao) = -4(2^4), P2M = 0. 

This Lemma gives the existence of a solution of a boundary value problem for 
the Hamiltonian system of equations with Hamiltonian T. It is this solution whose 
symmetric extension gives a symmetric periodic simultaneous binary collision orbit 
in the planar pairwise symmetric four-body equal mass problem. 

Theorem 5.2. Fix m= 1. For each E <0, there exists a time-reversible periodic 
regularized simultaneous binary collision orbit 

7(s) = {ui{s),U2{s),U3{s),U4{s),Vi{s),V2{s),V3{s),V4{s)) 

with period T > 0, angular momentum A = 0, and a symmetry group isomorphic 

to D4, for the Hamiltonian system of equations with Hamiltonian T on the level 
set r = such that distinct regularized simultaneous binary collisions occur at 
s = T/8, 3T/8, 5T/8, 7T/8. This periodic orbit corresponds to a symmetric periodic 
singular orbit 

{Xi{t), X2{t), X3{t), X4{t), UJi{t), u;2{t), u;3{t), u;4{t)) 

with energy E for the planar pairwise symmetric four-body equal mass problem 
where for all t, 

Xiit) = Xi{t), X3{t) = X2{t), \X2{t)\ < Xi{t), UJ4{t) = UJi{t), UJ^it) = U)2{t), 

with initial conditions 

xi{0) > 0, a;2(0) = 0, wi(0) = 0, ^2(0) > 0, 

and period 

.T/2 

R= / {ul{s)+ul{s)){ul{s)+ul{s)) ds, 
Jo 

where for t e [0,R], the only singularities are two distinct simultaneous binary 
collisions occurring att = R/A and t = 3-R/4. 
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Proof. Let Qi{(j), Q2i<^), -Pi(o'), P2{cr), cr G [0,(To], be the nonsingular solution of 
the Hamiltonian system of equations with Hamiltonian F on the level set F = 0, 
whose existence and properties are given in Lemma 5.1. Using the scaling a = s/jjb, 
set So = MCTo- By the linear symplectic transformation with multiplier we have 



21/4 ' -v"^ 2V4 ' 2^/^723l' 2^7^! 

Since «2 = — (v^— U4 — (a/2 — 1)^3, i'2 = — (v^— l)wi, and W4 = (-^2— l)w3, 
we have 

-(V2-1)Qi(s/m) . . (V2-1)Q2(s/m) 

-Vy2-lPi(s/M) . . Vy2-lP2(s//x) 
^^2(5) = ^ , Vi{s) = . 

Set 7(s) = {ui{s),U2{s),U3{s),Ui{s),v\{s),V2{s),V2,{s),Vi{s)), s e [0, So]- From 
£ = (192 _ 16^/2 - 8)/8 < and £; = (4 - 2\/2)£; and £; = 2E/{A - 2spif, we get 

16 ^-^<°- 

With this value of it follows that the value of F at 7(0) is 0. Set 

Ci = ^>0' C2 = ^^^f^<o, 

Pi = = < 0, P2 = ^^^ >0. 

/2 - 1 2 



With Qi(cr), Q2(o'), -Pi(o'), P2{p), cr € [O,cro], being a nonsingular solution of the 
Hamiltonian system of equations with Hamiltonian F, we have Q\{p^ + Q2{'^) 7^ 
for all a G [O,0-o]. From this it follows for all s G [0, so] that 

Ms = ul{s) - ulis) + ul{s) - ul{s) = {2-V2) [^^(s/m) + Ql{s/fi)] ^ 0, 

and 

Ms = 2ui(s)m2(s) - 2u3(s)u4(s) = -(2 - V2)[Qlis/fi) + Q^(s//z)] ^ 0. 

These imply that M| + M| 7^ and M| + M| ^ for all s G [0, sq]. Thus the 
function j{s) is a nonsingular solution of the Hamiltonian system of equations with 
Hamiltonian F on the level set F = that satisfies 

Wl(0)=Cl, W2(0)=C2 W3(0)=Cl, «4(0) = -C2, 

Wl(0)=pi, t;2(0) = P2, V3(0) = -pi, V4{0)=P2, 

wi(so)=0, W2(so)=0, ■U3(so)>0, W4(so) > 0, 

ui(so)<0, W2(so)>0, W3(so)=0, W4(so) = 0. 

By Lemma 4.1, the solution j{s) extends to a T = 8so periodic solution, call it 
7(s), with a D4 symmetry group generated by the symmetries Sp and Sg, and four 
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distinct regularized simultaneous binary collisions at s = Sq, Ssq, 5so, 7so, for which 

Wi(so) + "2(50) = 0, ul{so) + uliso) ^ 0, 

uf (3so) + wi(3so) 7^ 0, ul{3so) + uI{3sq) = 0, 

Wi(5so) + uli^so) = 0, ul{5so) + ^4(580) 7^ 0, 

uf (7so) + ui(7so) ^ 0, ul{7so) + ul{7so) = 0. 

Since Qi{a) + Qi{a) ^ for all a G [0,cto], it follows that (uf (s) + ui(s))(ui(s) + 
U4(s)) 7^ for s e [0,T] except s = so,3so,5so,7so. The regularizing change of 
time 

defines t an invertible differentiable function of s, i.e., t = 6{s) with ^(0) = and 

0'{s) = when s = {2k + l)so for fc G Z. The symmetry Sp satisfies Sf'^{s) = 
j{s + 2so) and —j{s) = Spj{s) = 7(s + 4so). The symmetry Sq satisfies 5^7(5) = 
7(2so — s), and so 7(s) has a time-reversing symmetry. The angular momentum of 
7(5) at s = is 

A = ^ [ - t^iU2 + V2U1 - V3U4, + U4U3] = P2C1 - P1C2 = 0. 
The extended solution 7(s) gives a singular symmetric solution 

Z{t) = {xi{t),X2{t),X3{t),X4{t),iOl{t),UJ2{t),u;3{t),U)4{t)), 

of the planar pairwise symmetric four-body equal mass problem. Under the Ansatz, 
the components of z{t) satisfies X4{t) = xi{t), X3{t) = X2{t), \x2{t)\ < xi{t), W4(<) = 
uJi{t), 0J3{t) = u)2{t), where 

^ ul{s)-ul{s)+ul{s)+ul{s) 

X2{t) = Ul{s)u2{s) +U3(s)u4(s), 

_ Vl{s)ui{s) - V2{s)u2{s) V3{s)u3{s) - V4{s)u4{s) 

'''^'>- 2iul{s)+ul{s)) + 2{uUs)+ulis)) ' 

_ Vl{s)u2{s) + V2is)ui{s) V3{s)u4{s) + V4{s)u3{s) 

2{ul{s)+ul{s)) + 2{ul{s)+ul{s)) ' 

for s = 9~^{t). The components of the extended solution 7(5) satisfy ^3(0) = Ui(0), 
W4(0) = -U2(0), ui(0)u2(0) < 0, |U2(0)| = (x/2- l)|ui(0)|, U3(0) = -?;i(0), Vii^) = 

V2{0), 

Vl(0)U2(0) + U2(0)ui(0) = P1C2 + P2C1 = ^1/4 > 0, 

and 



^^2(0)U2(0) - ui(0)ui(0) = P2C2 - PiCi = " 2^4^"^ - 

From Lemma 2.1, it follows that a;i(0) > 0, X2(0) = 0, uJi{0) = 0, and ^2(0) > 0. Set 
R = 9{T/2). Since 7(4so) = -7(0), it follows that xi{R) = xi{0), X2{R) = a;2(0), 
LJ\{R) = wi(0), and uj2{R} = ^^2(0). Thus the singular symmetric solution z{t) has 
period R. By the construction of the extension of 7(5) given in Lemma 4.1, there 
holds 

/■(fe+l)so rso 

/ (uf (s) + ui(s))(u^(s) + ul{s)) ds= (uf (s) + ul{s)){ul{s) + ul{s)) ds 

J ksQ Jo 
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for all A; = 1, . . . , 7. This implies that i?/4 = 0{{k + l)so) - d(kso) for all k = 
0, 1, . . . , 7. The first regularized simultaneous binary collision for 7(5) occurs at 
s = So, and this corresponds tot = 6{so) = R/A. The next regularized simultaneous 
binary collision for 7(5) occurs at s = 3so, and this corresponds to 

t = e(3so) = (^(3so) - e{2so)) + {0{2so) - e{so)) + e{so) = | + | + ^ = ^. 

Similarly, the regularized simultaneous binary collisions for 7(5) occurring at s = 
5so,7so correspond to t = 5i?/4, 7i?/4. Hence, for t E [0,R], the periodic solution 
z{t) has simultaneous binary collisions as its only singularities, and these occur at 
t = R/4, 3i?/4. 

For a fixed but arbitrary e > 0, the value of e~'^E is a fixed but arbitrary negative 
real number. By Lemma 3.1, the scaled extended solution 

7e(s) = (eui(es), eU2{es), eU3{es), eU4{es),Vi{es),V2{es), V3{es), v^ies)) 

is a periodic solution for the Hamiltonian system of equations with Hamiltonian T 

on the level set F = 0, having period e~^T and energy e~'^E < 0. The extended 
solution 7e(s) has four distinct regularized simultaneous binary collision occurring 
at s = {2k — l)e~^So for k = 1,2, 3, 4. The regularizing change of time now defines 
t as an invertible differentiable function of s by 

t = [ e\,4{ey) + uUey)){uliey) + ul{ey)) dy 
Jo 

e^uUy) + ul{y)){ul{y) + ul{y)) dy = e^e{es). 





The extended solution 7e(s) defines a singular symmetric solution 

z,{t) = {xl{t),xl{t),xl{t),xl{t),iol{t),u'2{t),u;l{t),ojl{t)) 

of the planar pairwise symmetric four-body equal mass problem, with xl{t) = xf (t), 
xl{t) = \xm\ < xlit), ujlit) = ul{t), = where 



x\{t) 



uf(es) — ulies) + ul{es) + ul{€s) 



2 

2^2(0 = e^['«i(es)M2(es) + U3{es)u4{es)], 

vi{es)ui{es) - V2{es)u2{es) ^ V3{es)u3{es) - V4{es)u4ies) 

LOi[t) — 



2t{u\{ts) ^ ul{^s)) 2e(ui(es)+u|(es)) 
_ V\{ts)u2{es) + V2{es)ui(ts) V3{es)u4{ts) + V4{es)u3{ts) 
2t{u\{ts) + ul{ts)) 2e{ul{ts) ^ u\{es)) 

for es = B-^{e-H). The period of z^{f) is e'R = €^0{T/2). The argument for the 
case of e = 1 now applies to e 7^ 1 with 7(5) replaced by 7e(s), Sq replaced by esg, 
z{t) replaced by Ze{t), T replaced by e~^T, E replaced by e^'^E, and R replaced by 
r'^R. This gives the existence of a symmetric periodic simultaneous binary collision 
orbit {t) of the planar pairwise symmetric four-body equal mass problem for every 
negative value of energy. □ 
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6. Numerical Estimates in the Equal Mass Case 

In the equal mass case, there is by Theorem 5.2 a time-reversible periodic orbit 
7(s) for the Hamiltonian system of equations with Hamiltonian T on the level set 
r = with period T. The components mi(s), U2{s), u^{s), U4(s), wi(s). f2(.s), ^3(3), 
V4{s) of 7(5) satisfy W2(s) = -(V2-l)ui(s), 144(5) = (a/2-1)w3(s), 112(5) = -(a/2- 
l)i;i(s), and ^4(5) = {V2 — l)i;3(s). The D4 symmetry group of 7(5) is generated 
by Sf"/{s) = 7(5 + T/4) and Sg"/{s) = 7(T/4 — 5). Under the linear symplcctic 
transformation with multiplier jj, this gives a periodic orbit Qi{a), Q2{o'), -F'i(c), 
^2(0") of the Hamiltonian system of equations with Hamiltonian F on the level set 
r = 0, which by Lemma 5.1 satisfies Qi{0) = 1, 02(0) = 1, Pi(0) = -1?, P2(0) = t9, 
Qi(cro) = 0, Q2{<Jo) > 0, Pi(c7o) = -4(2^/4), and P2(c7o) = for some CTo > and 
^? > 0. In [2], we numerically estimated 

(70 = 1.62047369909693, 1? = 2.57486992651942. 

The period of this periodic orbit is Sctq ~ 12.96378959 and its energy is i? sa 
—2.999682732. From the linear symplectic transformation with multiplier and 
the relations among the components of 7(5), we have for the values of components 
of 7(0) the exact 

ui(0) = 1x3(0) = 2-V4, -U2{0) = 1x4(0) = (\/2 - 1)2-1/4, 

and the estimates 

i;i(0) = -i;3(0) = , ss -2.000382939, 

2V^/2- 1 

V2{0) = V4{0) = -^-^ « 0.8285857433. 

Since a = s/fx, the period of 7(5) is T « 8.469003682. Since E = E/{2 - \/2), the 
value of the energy for 7(5) is ^ « -5.120778733. 

Figure 2 illustrates the graphs of the components of the scaled periodic orbit 7e(5) 
for an e > 2. Readily observable in these graphs arc the symmetries Spjeis) = 
7e(5 + T'j/4) and Sg%{s) = 7e(Tj/4 — 5). Figure 1 illustrates the curves in the 
physical plane that the four equal masses follow in the symmetric simultaneously 
binary collision orbit z{t) = {xi{t), X2{t), 2:3 (t), X4{t), uj{t), 0^2 (i), 1^3 (i), <^4{t)) 
of the planar pairwise symmetry four-body problem, corresponding to % (5) for an 
e > 0. The initial conditions for z{t) are 

a;i(0) = a;4(0) = 1, a;2(0) = ^3(0) = 0, 

a;i(0) = a;4(0) = 0, ^2(0) = ^3(0) w 1.287434964. 

The value of e = l/\/2 — \/2 here for the scaling is determined by the relation 
Xi{t) — e^(itj(es) — U2{es)) coming from the regularizing transformation applied 
to the scaled periodic solution 7e(s), together with the initial condition a;i(0) = 1, 
where t = corresponds to s = 0. The value of the Hamiltonian H along z{t) is 
e-'^E w -2.999682732. 

Some analytic evidence for the linear stability of 7(5) is provided by an inves- 
tigation of the linear stability of the periodic orbit Qi{<t), Q2{cr), Pi(f)) P2(o') 
of the Hamiltonian system of equations with Hamiltonian F, corresponding to 
71(5) = 7(5). In [2], we showed that the periodic orbit Qi{(t), (52(f), Pi(o'), P2(c) 
is linearly stable for the Hamiltonian system of equations with Hamiltonian F. We 
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Figure 2. The ui, it2, U3, it4, vi, V2, v^, V4 coordinates of the 
symmetric periodic simultaneous binary coUision orbit 7e(s) for 
m = 1 and an e > 2. 




Figure 3. A long-term numerical integration of the symmetric 
periodic simultaneous binary collision orbit 7e(s) for m = 1 and 
an e > 2. 



did this by applying the symmetry reduction technique of Roberts [16] to compute 
its two nontrivial characteristic multipliers to be a complex conjugate pair lying on 
the unit circle close to but not equal to —1. (Of course, the trivial characteristic 
multiplier for this periodic orbit is 1 whose algebraic multiplicity is two). This 
means that 7(5) is linearly stable among all solutions of the Hamiltonian system 
of equations with Hamiltonian F on the level set F = whose components satisfy 
U2 ~ -(\/2 - l)ui, U4 — {V2 - l)u3, V2 — -(\/2 — and V4 — {V2 — 1)1)3. 
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Numerical evidence for the linear stability of the periodic orbit 7(5) is provided 
by a long-term numerical integration of the Hamiltonian system of equations with 
Hamiltonian F. We used the Runge-Kutta-Fehlberg algorithm, with a time-step of 
0.001, for the numerical integration. By Theorem 3.3, the linear stability of 7e(s) 
is the same for all e > 0. Figure 3 illustrates the components of 7e(s), for the same 
e > 2 as in Figure 2, over the time interval [0,100], where 100 is approximately 
26re. We also performed a longer term numerical integration of 7e(s) for the value 
of e that satisfies — I-k (see next Section). From these numerical integrations it 
appears that 7e(s) is linearly stable for some value of e, and hence for all e > 0. 

Further numerical evidence of the linear stability of the periodic orbit 7(5) is 
provided by a long-term numerical integration of the orbit of four equal masses 
starting at the initial conditions determined by 2;(0). We did this, not for the 
Hamiltonian system of equations with Hamiltonian iJ, but for the three-dimensional 
four-body equal mass problem. We carried out the numerical integration using 
Euler's method combined with a subroutine that gives an elastic bounce at binary 
collisions. None of the symmetry properties of the singular periodic simultaneous 
binary collision orbit nor the pairwise symmetry of the four masses, were coded 
into the algorithm. Because the motion of the four masses is planar, the numerical 
integration produces an orbit of the planar four-body equal mass problem. Any 
instability in this planar orbit should be amplified by Euler's method and become 
apparent in a timely manner. The numerical integration from f = to i « 100000, 
with a time-step of 0.001, shows the orbit remains bounded, with the four masses 
tracing out planar curves close to those illustrated in Figure 1, while retaining its 
regular pattern of alternating simultaneous binary collisions. 

7. Numerical Estimates in the Unequal Mass Case 

We numerically continue to < to < 1 the time-reversible periodic regularized 
simultaneous binary collision orbit for the Hamiltonian system of equations with 
Hamiltonian F on the level set F = 0, whose analytic existence is given by Theorem 
5.2 for TO = 1. We also investigate through long-term integrations the values of 
< TO < 1 for which the continued timc-revcrsiblc periodic regularized simultaneous 
binary collision orbits are linearly stable. For to = 1, we assume by Lemma 3.1 that 
the time-reversible periodic regularized simultaneous binary collision orbit 7(5; 1) 
with the D4 symmetry has period T = 2tt and energy E w —2.818584789. Recall 
that the D4 symmetry group of 7(5; 1) is generated by Sf^{s; 1) = 7(s-|-7r/2; 1) and 
'Salis': 1) = 7(7r/2 — s; 1). We assume without loss of generality by Lemma 3.1 and 
Lemma 4.1, that the continued time-reversible periodic regularized simultaneous 
binary collision orbit "i{s; m) for < m < 1, has period T — 27r, energy E{m), and 
a D4 symmetry group generated by 5^7(5; to) = 7(5 -|- Tr/2;m) and S'g7(s;to) = 
7(71/2 — s; to). We assume that E{m) is a continuous function of < to < 1, with 
E{1) = E. We shift the regularized time variable s to s -|- 7r/4, so that for 

7(s; to) =(tti(s; to), 1*2(5; to), ^3(5; to), ^4(5; to), 
vi{s; m),V2{s; to), ^3(5; to), V4{s; to)), 

the value s = now corresponds to the first simultaneous binary collision, i.e., 
ui(0;to) = 0, M2(0;to) = 0, W3(0;to) ^ 0, M4(0;to) ^ 0, vi{0;m) ^ 0, V2{0;m) ^ 0, 
W3(0;m) = O, and V4{0,m) = 0. We approximate the continued time- reversible 
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periodic orbits 7(5; m), < m < 1, by the trigonometric polynomials, 



n n 



ui{s;m) = tti sin((2i — l)s), U2{s; m) = — 6j sin((2i — l)s) 




U3(s; m) = ^1(8 — 7r/2; m), ^4(5; m) = U2{s + 7r/2; m), 



n n 



vi{s; m) = Cj cos((2« — l)s), 1^2(5; m) = — cos((2i — l)s) 




^3(5; m) = Wi(s — 7r/2; m), W4(s;m) = V2{s + 7r/2; m), 

for a positive integer n and constants ai, bi, Ci, di, i = 1, . . . ,n, that are assumed to 
be continuous functions of m. The presence of odd positive integer frequencies 2i — l 
in these trigonometric polynomials is to ensure that the period functions defined 
by them have the D4 symmetry group generated by Sp and Sq- So in particular, 
the periodic orbits in the continuation are time- reversible for all < m < 1. An 
numerical estimate of the periodic solution j{s; m) is found through the variational 
approach of minimizing the functional 



over the space M of coefficients a^, 6^, c^, d^, z = 1, . . . , n, for an appropriate choice 
of n. The use of trigonometric polynomials for numerically approximating periodic 
solutions is a classic approach (see, for example, Simo [20]). 

The mimcrical algorithm for finding a trigonometric polynomial that approxi- 
mates the periodic solution 7(3; to), < to < 1, proceeds in two steps. First, we 
consider a guess for the set of values for a,, bi, Ci, di, i = 1, . . . , n, as well as a guess 
for E{m). Starting with a reasonably low number of terms, n = 5, we let a numer- 
ical minimization algorithm (in this case, MATLAB's fminunc) find the minimum 
of L near the starting guess. Then we add an additional non-zero term to each 
of the trigonometric polynomials, and the minimizing solution from the previous 
iteration is used as a starting guess for the next iteration. This process continues 
until we reach n = 10. Second, since the Hamiltonian equations with Hamiltonian 
r requires a specified value of E{m) to compute L, we need to make certain that 
we are getting a good estimate of E{m). We evaluate F at the point 7(7r/4;TO), 
away from simultaneous binary collisions. If this value is not sufficiently close to 
(within about 5 x 10"^^° of 0), we adjust E{m) using the bisection method in a 
small interval about the initial guess of E{m) until r(7(7r/4; to)) is sufficiently close 
to 0, re-minimizing the trigonometric polynomial approximation of 7(5; m) for each 
new choice of Eim). 

This numerical method only works well if we have a good initial guesses for Oj, 
bi, Ci, di, i = 1, . . . ,5, and E{in) for some value of m. This we have when to = 1. 
We use our estimate of 7(5; 1) and E(l) to provide the initial guesses for a,, 6,, 
Ci, di, i = 1, . . . ,5, and E{m) for to = 0.99. The mimcrical algorithm produces a 
trigonometric polynomial approximation of 7(5; 0.99) and an estimate of i?(0.99), 
which then provide the initial guesses for Oj, bi, Ci, di, i = 1, . . . ,5, and E{m) for 
TO = 0.98. We continue decreasing to by 0.01 and using the numerical algorithm 
until we reach to = 0.01. In Figure 4, we plot the graphs of the numerical estimates 
of U3(0;to), 'U4(0;to), wi(0;to), and V2{0;m). Notice that Vi{0;m) + V2{Q;m) w 
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Figure 4. The estimated values of 7/3(0; m) (blue), U4(0;m) 
(green), wi(0;m) (red), and V2{0;m) (light blue) for the nonzero 
components of 7(6; to) over < m < 1. 
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Figure 5. The energy E{m) of the periodic orbits 7(5; m) for 

< 771 < 1. 

32to^/(to + 1), as is expected from the regularization of the simultaneous binary 
collisions. In Figure 5, we plot the graph of the numerical estimate of E{m). In 
Figure 6, we graph of value of r(7(7r/4; 777)) over < m < 1. 
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Figure 6. The value of f (7(71/4; to)) for the optimized trigono- 
metric polynomial approximation of 7(5; to) over < m < 1. 

We use the optimized trigonometric approximations of 7(5; to), < to < 1, to 
estimate 7(— 7r/4;TO), which are the conditions corresponding to s = before the 
shift of s. For each < m < 1, the conditions 7(— 7r/4; to) satisfy Lemma 2.1, and 
so they corresponds to the initial conditions of interest as given in the Introduction. 

We also use the optimized trigonometric approximation of 7(5; to) to estimate 
7(0; to). Using the estimated value of £'(to), for values of to between 0.01 and 1 
separated by 0.01 increments, we integrated the Hamiltonian system of equations 
with Hamiltonian T using the Runge-Kutta-Fehlberg algorithm, with a time-step 
of 0.001, starting at 7(0; to). By the time s reached 407r, i.e., twenty periods, the 
periodic orbits were easily recognized as unstable for m < 0.53. On the other 
hand, with s reaching values much larger than 407r, the periodic orbits appear to 
be linearly stable for to > 0.54. 

Following these observations, we refine the numerical continuation to approxi- 
mate the periodic orbits 7(3; to) and the energy E{m) for to = 0.531, 0.532, . . . , 
0.539 using the numerical algorithm described previously. Then the Hamiltonian 
systems of equations with Hamiltonian F were integrated by the Runge-Kutta- 
Fehlberg algorithm, with a time-step of 0.001, for each of the initial conditions 
7(0; to) with m = 0.531, 0.532,..., 0.539, until the norm of the orbit, ||7(s;TO)|j, 
exceeded 100 for some s > 0. By the time s reached IOOOOtt, i.e, 5000 periods, the 
only orbits whose norms were still below 100 were those with 0.535 < to < 0.539. 
The norm of the orbit for to = 0.535 exceeded 100 after s reached 156887r, i.e., 7844 
periods. The integrations for the orbits with masses 0.536 < m < 0.539 continued 
until the value of s was approximately 2018647r, i.e., 100,932 periods, at which 
point the integrations were halted. 

The graphs of the numerical integration of the components of 7(3; 0.536) are 
illustrated in Figure 7, where s = on the horizontal axis corresponds to s ~ 
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Figure 7. Graphs of the components of 7(5; m) when m = 0.536, 
where s = on the horizontal axis corresponds to s 2018647r. 




Figure 8. The motion of the four masses zt{xi{t), X2{t)) and 
±{x:}{t),Xi{t)) when m — 0.536 for large t, corresponding to the 
numerical integration of 7(5; 0.536) with s sa 20 1 8647r. 

2018647r. These graphs clearly show that the numerical integration of 7(5; 0.536) 
has ceased to have a period of 27r. Illustrated in Figure 8 is the motion of the 
four pairwise symmetric masses corresponding to the numerical integration of the 
components of 7(5; 0.536) as shown in Figure 7. The rotational drift from the initial 
conditions and the irregularity of this motion indicates that there is instability near 
the periodic orbit 7(5; 0.536). 

The graphs of the numerical integration of the components of 7(3; 0.537) are illus- 
trated in Figure 9, where s = on the horizontal axis corresponds to s « 2018647r. 
From these graphs it appears that the numerical integration has retained the peri- 
odicity of 27r. Illustrated in Figure 10 is the motion of the four pairwise symmetric 
masses corresponding to the numerical integration of the components of 7(5; 0.537) 
as shown in Figure 9. The rotational drift from the initial conditions indicates 
that there is an instability near the periodic orbit 7(5; 0.537). The regularity of the 
motion in Figure 10 indicates that m = 0.537 is close to the value of m at which 
the 7(3; m) is spectrally stable but not linearly stable. 

The graphs of the numerical integrations of the components of 7(3; 0.538) are 
illustrated in Figure 11, where s = on the horizontal axis corresponds to s sa 
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Figure 10. The motion of the four masses zt{xi{t),X2{t)) and 
±{x3{t),X4{t)) when m — 0.537 for large t, corresponding to the 
numerical integration of 7(5; 0.537) with s « 2018647r. 



2018647r. From these graphs it appears that the numerical integration of 7(5; 0.538) 
has retained the periodicity of 2tt. Illustrated in Figure 12 is the motion of the 
four pairwise symmetric masses corresponding to the numerical integration of the 
components of 7(3; 0.538) as shown in Figure 11. The lack of rotation and the 
presence of a thin annulus-like region traced out by the motion of four pairwise 
symmetric masses indicates that 7(5; 0.538) is linearly stable. 

The graphs of the numerical integrations of the components of 7(5; 0.539) are 
illustrated in Figure 13, where s = on the horizontal axis corresponds to s ~ 
2018647r. From these graphs it appears that the numerical integration of 7(5; 0.539) 
has retained the periodicity of 2tt. Illustrated in Figure 1 is the motion of the 
four pairwise symmetric masses corresponding to the numerical integrations of the 
components of 7(3; 0.539) as shown in Figure 13. The lack of rotation and the 
regularity of the motion indicates that 7(3; 0.539) is linearly stable. 

This numerical investigation shows that 7(5; 1) continues as a one-parameter 
family 7(5; m), < to < 1, of symmetric periodic simultaneous binary collision 
orbits, all with period 27r and varying energy E(m). It further shows that 7(5; to) 
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Figure 11. Graphs of the components of 7(5; m) when m = 0.538, 
where s = on the horizontal axis corresponds to s 2018647r. 




Figure 12. The motion of the four masses zt{xi{t),X2{t)) and 
zL{x^{t),X4{t)) when m = 0.538 for large t, corresponding to the 
numerical integration of 7(5; 0.538) with s 2018647r. 




Figure 13. Graphs of the components of 7(5; m) when m — 0.539, 
where s = on the horizontal axis corresponds to s sa 2018647r. 



is linearly stable when 0.538 < m < 1 and is unstable when < m < 0.537, with 
spectral stability occurring without linearly stability for a value of m close to 0.537. 
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EXISTENCE AND STABILITY OF SYMMETRIC PERIODIC 

SIMULTANEOUS BINARY COLLISION ORBITS IN THE 
PLANAR PAIRWISE SYMMETRIC FOUR-BODY PROBLEM 



LENNARD F. BARKER, TIANCHENG OUYANG, DUOKUI YAN, AND SKYLER SIMMONS 

Abstract. Wc extend our previous analytic existence of a symmetric periodic 
simultaneous binary collision orbit in a regularized fully symmetric equal mass 
four-body problem to the analytic existence of a symmetric periodic simulta- 
neous binary collision orbit in a regularized planar pairwise symmetric equal 
mass four-body problem. We then use a continuation method to numerically 
find symmetric periodic simultaneous binary collision orbits in a regularized 
planar pairwise symmetric 1, m, 1, m four-body problem for m between and 
1. Numerical estimates of the the characteristic multipliers show that these 
periodic orbits axe linearly stability when 0.54 < m < 1, and are linearly 
unstable when < m < 0.53. 



1. Introduction 

In the A/'-body problem, linearly stable periodic orbits may trap around them- 
selves bounded, non-chaotic motion of the N masses [5]. Some of the known ex- 
amples of linearly stable periodic orbits in the three-body problem are the elliptic 
Lagrangian triangular periodic orbits for certain values of eccentricity and the three 
masses [12], [18], and the Montgomery-Chenciner figure-eight periodic orbit for equal 
masses [4], [14], [19], [8], [9]. Other examples of linearly stable periodic orbits in the 
three or four-body problem involve binary collisions (BC) and/or simultaneous bi- 
nary collisions (SBCs). The rcgularization of these kinds of singularities is achieved 
by a generalized Levi-Civita type transformation and an appropriate scaling of time, 
as adapted from Aarseth and Zare [1] to the particular problem (also see [10], [16]). 

Schubart [26] numerically discovered a singular symmetric periodic orbit in the 
collinear three-body equal mass problem. In this orbit, the inner mass alternates 
between binary collisions with the two outer masses. Henon [6] extended Schubart's 
numerical investigations to the case of unequal masses. Only recently did Vonttirclli 
[29] and Moeckcl [13] prove the analytic existence of the Schubart orbit when the 
outer masses arc equal and the inner mass is arbitrary. The linear stability of the 
Schubart orbit was determined numerically by Hietarinta and Mikkola [7] revealing 
that linear stability occurs only for certain choices of the three masses. Numerically, 
non-Schubart-like linearly stable periodic orbits in the collinear three-body problem 
were found by Saito and Tanikawa for certain choices of the masses [20], [21], [22]. 

Sweatman [27], [28] and Sekiguchi and Tanikawa [24] numerically found a sym- 
metric Schubart-like orbit in the symmetric collinear four-body problem with masses 
1, m, m, and 1. This Schubart-like periodic orbit alternates between SBCs of the 
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70H33. 
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two outer pairs of masses and binary collisions of the inner two masses. Ouyang 
and Yan [17] proved analytically the existence and symmetry of this orbit. In the 
regularized setting, this periodic orbit has a symmetry group isomorphic to D2, 
of which both of the generators are time-reversing symmetries. (Here the dihedral 
group Dfc is the group of symmetries of the regular fc-gon.) Swcatman [28] numer- 
ically showed that the Schubart-like orbit is linearly stable when < m < 2.83 or 
m > 35.4, and is otherwise unstable. This linear stability was confirmed [3] using 
Robert's symmetry reduction technique [19]. 

Ouyang, Yan, and Simmons [15] numerically found and analytically proved the 
existence and symmetries of a singular symmetric periodic orbit in the fully sym- 
metric planar four-body problem with equal masses. (In the fully symmetric planar 
four-body equal mass problem, the position of one mass determines the positions of 
the other three masses.) In this orbit, the four masses alternate between different 
SBC's. In the regularized setting, this periodic orbit has a symmetry group isomor- 
phic to D4, of which one of the two generators is a time-reversing symmetry, while 
the other generator is a time-preserving symmetry. By using Roberts' symmetry 
reduction method, we showed [3] that this symmetric periodic simultaneous binary 
collision orbit is linearly stable. 

In this paper we extend the analytic existence of a symmetric periodic SBC orbit 
in the fully symmetric planar four-body equal mass problem [15] to the analytic 
existence of a symmetric periodic SBC orbit in the planar pairwise symmetric four- 
body problem, or PPS4BP for short. The positions of the four pairwise symmetric 
bodies in the plane are (xi,X2), {x^jX^), {—Xi,—X2), and (— ^3, — .T4), where the 
corresponding masses are 1, m, 1, m with < m < 1. With t as the time variable 
and ■ = d/dt, the momenta for the four masses are (a;i,W2) = 2(xi,X2), (w3,W4) = 
2m(i3,a;4), —{ijJ\,uj2), and —{uj^^uji). The Hamiltonian for the PPS4BP is if = 
K — U, where 



and 

U 



1 2m 



2^x1 + x\ sj{x-i - Xif {XA - X2Y 

2m m^ 
+ .„ , -I- 



^{xx xzf {X2 x^f 2^x1 x\ ■ 
The angular momentum for the PPS4BP is 

A = X\ljJ2 — X2OJ1 + X3LO4 — XiU)^. 

The center of mass is fixed at the origin, and the linear momentum is zero. With 



J = 



I 
-I 



for I the 4x4 identity matrix, the vector field for the PPS4BP is JVH, i.e., the 
Hamiltonian system of equations with Hamiltonian H are Xi = dH/dui, tbi = 
—dH/dxi, i = 1,2,3,4. The PPS4BP presented here is the Caledonian symmetric 
four-body problem [23] with non-collinear initial positions. 

The initial conditions for the orbits of interest has the first body of mass 1 located 
on the positive horizontal axis with its momentum perpendicular to the horizontal 
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Figure 1. The symmetric periodic SBC orbit in the PPS4BP for 
m = 1 (left) and m = 0.539 (right). The two red curves are those 
traced out by ±{xi{t),X2{t)), and the two blue curves are those 
traced out by ±{x3{t), X4{t)). 



axis, and the first body of mass m located on the positive vertical axis with its 
momentum perpendicular to the vertical axis. Specifically, at t = we have 

xi > 0, X2 = 0, ^3 = 0, X4 > 0, with X4 < xi, 
oji =0, UJ2 > 0, W3 > 0, U4 — 0, with u!2 < a;3, 

at which H is defined. The first objective is to find, for < m < 1, values of 
xi, X4,u!2,i^3 at t = such that (i) x^ — xi = and X4 — X2 — with + X2 7^ 
at some t — to > 0, (ii) xi + x^ — and X2 + X4 — with xf + X2 ^ al some 
t = ti > ioi (fii) the orbit extends to a symmetric periodic orbit, and (iv) the 
periodic orbit avoids all the other kinds of collisions. Such an orbit experiences a 
SBC in the first and third quadrant at t = tQ, and then another SBC in the second 
and fourth quadrants at t = ti, before returning to its initial conditions at some 
t = t2 > ti. The presence of collisions along the orbit necessarily imposes zero 
angular momentum on the orbit, thus requiring that X1UJ2 ~ X4UJ3 = at t = 0. 
Examples of these symmetric periodic SBC orbits in the PPS4BP with masses 1, m, 
1, m are illustrated in Figure 1 for m — I and m — 0.539. The second objective is 
to numerically investigate the linear stability of the symmetric periodic SBC orbits 
as m varies over interval (0, 1]. 

The regularization of the SBCs, as described by (i) and (ii) above, in the Hamil- 
tonian system of equations with Hamiltonian H plays a key role in achieving the 
two objectives. Section 2 details this regularization which consists of two canonical 
transformations followed by a scaling of time t — 9{s) with s as the regularizing time 
variable, producing a new Hamiltonian F for the PPS4BP in extended phase space. 
Section 3 describes a scaling of orbits of the Hamiltonian system of equations with 
Hamiltonian F which shows that any such periodic solution always belongs to a 
one-parameter family of periodic solutions for which the linear stability is the same 
for all periodic solutions in the family. Section 4 describes the symmetries of the 
Hamiltonian F which are used to construct periodic solutions with a D4 symmetry 
group generated by a time-reversing symmetry and a time-preserving symmetry. 
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In Sections 5 and 6, we prove the analytic existence of a periodic SBC orbit 7(s), 
with a D4 symmetry group, for the Hamiltonian system of equations with Hamil- 
tonian F with m = 1, and numcricaUy investigate its stabihty. The proof extends 
the analytic existence of a symmetric periodic SBC orbit in the fully symmetric 
planar four-body equal mass problem, as found in [15], to the PPS4BP with equal 
masses. Our numerical estimates of the characteristic multipliers show that this 
periodic orbit is linearly stable. 

In Section 7, we numerically continue the linearly stable symmetric periodic SBC 
orbit 7(5) from m = 1 to a symmetric periodic SBC orbit 7(5; m) for m < 1, and 
then investigate the linear stability of 7(s; m) as m varies in the interval (0, 1]. We 
use trigonometric polynomials as approximations of the periodic orbits (cf. [25]). 
The numerical algorithm for continuation starts with a trigonometric polynomial 
approximation of 7(s; 1) = 7(3) that is used as a initial guess for a trigonometric 
polynomial approximation of 7(5; 0.99), where the coefficients of the trigonometric 
polynomial are optimized through a variational approach. This process is repeated, 
using the optimized approximation of 7(s;0.99) as the initial guess for 7(5; 0.98), 
etc., until an optimized approximation of 7(5; 0.01) is obtained. Numerical esti- 
mates of the characteristic multipliers of 7(s; m) show it is linearly stable when 
0.54 < TO < 1 and unstable when < m < 0.53. We have analyzed the linear 
stability of these orbits for < to < 1, especially for 0.53 < m < 0.54 with incre- 
ments of 0.001, using Roberts' symmetry reduction method, the results of which 
will appear in a subsequent paper [2] . 

2. Regularization 

We adapt the regularization of Aarseth and Zare [1] to the PPS4BP to regularize 
SBCs as described in the first objective. This regularization differs from the one 
used in the Caledonian symmetric four-body problem [23] in that we only regularize 

the SBCs as described in the first objective of the Introduction. 
The first canonical transformation in our regularization is 

determined by the generating function 

Fi{xi,X2,X3, Xa, hi,h2, hz, hi) = hi{xi-X3) + h2{x2-X4) + h3{Xi+X3) + h4,{x2+X4,). 

So the first canonical transformation is determined by 

dFi dFi 
(1) ^^^~dx ' ^'^aoT' ^^l'^'^'"*- 

The new Hamiltonian is H = K — U, where 

- _ {hi + hs)^ + {h2 + /t4)' , (/13 - /ll)' + (^4 - h2)^ 

4 + 4to 

and 

1 2to 2to to^ 

U= ,„ , .„ + , „ . . „ + 



V(5i +53)2 + (52 + 54)' \^g!Tgl ^/ai + gi V(5i - 53)2 + (92 - 54)' ' 

The second canonical transformation in our regularization, 

(5l>fl'2,53,54,/ll,/l2,/l3,/l4) {u\,U2.,Uz,U4,V\,V2,V3,V^ 
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is determined by the generating function 



F2{hi,h2,h3,h4,Ui,U2,U3,U4) = - '^hjfj{ui,U2,Uz,Ui) 

where 

fx=ul-ul, f2 = 2uiU2, f3=ul-ul, f4 = 2usU4. 

So the second canonical transformation is determined by 
The new Hamiltonian is H = K ~ U, where 



dF2 



Id \ TO 



1 / J_\ (f3M3 - V4U4,){viUi - V2U2) + iv3U4 + 7747^3) (^1^2 + ^2^l) 



TO / 



(u2 + «2)(y2^y2) 



and 



2to 



2to 



TO^ 

+ 



+ + {2uiU2 - 2U3U4)2 



We introduce a new time variable s by the regularizing change of time 

^^={ul + ul){ul + ul). 

To simphiy notation, we set 

Ml = viui — V2U2, 

M3 = V^U^ - V4U4, 

M5 = u\ — u\ + u\ — 



M7 = u\ — u\ — u 



J- 4, 



M2 = V1U2 + V2U1, 

M4 = V^U4 + 'i-'4?i3, 
Affj = 2uiU2 + 2'U,3M4, 

Mg = 2uiU2 — 2M3U4. 



The Hamiltonian in the extended phase space with coordinates ui, U2, Us, U4, E, 
Vl, V2, V3, V4, t is 



dt 
ds 



H-E) 



1 

16 
+ 



1 + 



1 



1 - 



m 



{vl + vDiul + ul) + (vl + vl){ul + ul) 



(M3M1 + M4M2) 



li + Mi 

^2/„,2 I „,2\U,2 



2m + M2 + M3 + U4) 



m^{u{ +ui){u% +ui) 



E{u\ + ul){ul+ul). 



VM^ + Mi 

With ' = d/ds, the Hamiltonian system of equations with Hamiltonian T is 



(3) 



u'-^ v'--^ i-1234 
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along with the auxihary equations, 

(4) ^'=f =0' t' = -^ = {ul+ul){ul+ul). 

On the level set T = 0, the value of the Hamiltonian H (i.e., the energy) along 
solutions of the Hamiltonian system of equations with Hamiltonian T is E. Indepen- 
dent of the values of E and F, the angular momcntimi A = xiui2—X2(^i+X3U)4—XiL03 
in the coordinates ui,U2,U3,U4,vi,V2,V3, va simplifies to 

(5) ^ ^ ^ I- ~ ^^^^ ^^"^ ~ ^^^^-1 ■ 

On the level set F = 0, two simultaneous binary collisions in the PPS4BP have 
been regularized in the Hamiltonian system of equations with Hamiltonian F. The 
simultaneous binary collision — a;i = and .T4 — .X2 = with x\ + x\ 7^ 
corresponds to ul + u\ = Q with ul + u\ ^ 0. These imply that M| + M| = 
A{xl +xl)i^Q and + = 4(a;? + xl) 7^ 0. Prom f = it follows that 

(6) vl+vl = (Z2m^)/{m, + l). 

Similarly, the simultaneous binary collision X3 + a;i = and X4 + X2 = with 
x\ + x\ ^ corresponds to ul + ul = with ul+u^ 0, and hence that M^ + Mq = 
4(a;f + xl) 7^ 0, M| + M| = 4(a;? + x^) 0, and, from f = 0, that 

?;| + t;| = (32m2)/(m+l). 

On the level set F = 0, the other singularities of the PPS4BP have not been 
regularized in the Hamiltonian system of equations with Hamiltonian F. The binary 
collision x\ = 0, X2 = with x\ ^ x\ ^ ^ corresponds to M| + M| = and 
M| + M| 7^ with u\^u\^^ and uj + 0. The binary collision X3 = 0, 

X4 = with xl+ xl corresponds to M| + M| 7^ and M| + M| = with 
uf + ul and Wg + 7i| 7^ 0. Because of the pairwise symmetry, there are no triple 
collisions. Total collapse xi = 0, X2 = 0, X3 = 0, 2:4 = corresponds to ui = 0, 
U2 = 0, M3 = 0, and U4 = 0. A solution of the Hamiltonian system of equations with 
Hamiltonian F is called nonsingular if it avoids the unregularized binary collisions 
and total collapse singularities. 

We establish next the correspondence between the original coordinates and the 
regularized coordinates for the initial conditions given in the Introduction, the proof 
of which is straight-forward. 

Lemma 2.1. The conditions {att = 0) 

Xl > 0, X2 = 0, X3 = 0, X4 > 0, with X4 < Xl, 
oji =0, LO2 > 0, W3 > 0, W4 = 0, with u)2 < W3, 

correspond to the conditions {at s = 0) 

U3 = ±ui, U4 = with U1U2 < 0, \u2\ < (\/2 — l)|ui|, 
vs = TVi, V4 = ±V2, with < V1U2 + V2U1 < V2U2 — viui. 
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3. A Scaling of Periodic Orbits and Linear Stability 

A certain scaling of solutions of the Hamiltonian system of equations with Hamil- 
tonian T produces more solutions. When applied to a periodic solution, this scaling 
leads to a one-parameter family of periodic solutions. The proof of the following 
result is a straight-forward verification. 

Lemma 3.1. If"f(s) = {ui{s),U2{s),U3{s),Ui{s),vi{s),V2{s),V3{s),Vi{s)) is ape- 
riodic solution of the Hamiltonian system of equations with Hamiltonian T on the 
level setr = with period T and energy E, then for every e > 0, the function 

7e(s) = (eui(es), eU2{es), eus{es), eUi{es),v\{es),V2{es), Vs{es), Vi{es)) 

is a periodic solution of the Hamiltonian system of equations with Hamiltonian T 
on the level set F = with period = e^^T and energy = e~'^E. 

The linear stability of a periodic orbit of the Hamiltonian system of equations 
with Hamiltonian T is determined by the linearization of the equations along the 
periodic orbit. By Lemma 3.1. a periodic orbit 7(5) on the level set F = with 
period T and energy E embeds into a one-parameter family 7e(s) of periodic orbits 
on the level set F = with period and energy E^. The linearization of the 
Hamiltonian system of equations (3) with Hamiltonian F along the periodic orbit 
7e(s) is 

X' = JV^f (7,(s))X 

where V^F is the matrix of second-order partials of F. Let X^{s) be the solution 
of the linearization of the equations along 7e(s) that satisfies Xt (0) = / (the 8x8 
identity matrix). The monodromy matrix for 7e(s) is X^iTf,), and the eigenvalues 
of Xf_{Tf) are the characteristic multipliers of 7e(s). A characteristic multiplier 
A of 7e(s) is defective if its geometric multiplicity is smaller than its algebraic 
multiplicity, i.e., its generahzed eigenspace Uj>iker(Xe(Te) — \iy is not the same 
as its eigenspace ]ier{Xf {T^) — XT). Proving the following result is routine (see [11]). 

Lemma 3.2. If"f{s) is a periodic orbit of the Hamiltonian system of equations with 
Hamiltonian F on the level set F = 0, then for each e > 0, the periodic orbit 7e(s) 
has 1 as a defective characteristic multiplier with algebraic multiplicity at least two. 

For each e > 0, the periodic orbit "fe{s) is spectrally stable if all of its character- 
istic multipliers have modulus one. By Lemma 3.2, the periodic orbit 7e(s) has 1 
as a defective characteristic multiplier with algebraic multiplicity at least two, and 
so the monodromy matrix Xf{T^) it not semisimple. However, the two-dimensional 
subspace 

t^i=Span (7;(0),(^|:7.)(0)) 

is Xe(Te)-invariant. The periodic orbit 7e(s) is said to be linearly stable if it is 
spectrally stable and there exists a 6-dimensional Xg(Te)-invariant subspace U2 
such that Ui + U2 = M.^ and X^{Tf) restricted to U2 is semisimple. A proof of the 
following is routine. 

Theorem 3.3. Suppose 7(5) is a periodic orbit of the Hamiltonian system of equa- 
tions with Hamiltonian T on the level set F = 0. Then 7e(s) is spectrally (linearly) 
stable for some e> if and only if 'yds) is spectrally (linearly) stable for all e > 0. 
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4. Symmetries 



The Hamiltonian system of equations with Hamiltonian f has a group of sym- 
metries isomorphic to the dihedral group D4 = {a,b : = = (a6)^ = e). With 



-1 




G 



1 
1 



define the matrices 






F 








-F 




















F 








-F 






-G 














G 














G 














-G 



Sg 



These matrices satisfy Sp = —I, Sp = I, Sq = I, and {SfSg}'^ = I- Fixing the 
value of E, these matrices satisfy T o Sp = and F o S'g = F, and so Sp and Sq 
are the generators of the Z)4-symmetry group for F. If 

7(5) = {ui{s),U2{s),U3{s),U4{s),Vl{s),V2is),V3{s),V4{s)) 

is a solution of the Hamiltonian system of equations with Hamiltonian F, then 
5f7(s), 51.7(5), and 5g7(— s) are also solutions of the Hamiltonian system of 
equations with Hamiltonian F. This means that Sp is a time-preserving symmetry 
and that Sq is a time-reversing symmetry. The proof of the following is routine. 

Lemma 4.1. If for some sq > there is a nonsingular solution 7(s), s G [0, sq], 
of the Hamiltonian system of equations with Hamiltonian F such that for constants 
Ci 7^ 0; C2 7^ 0, pi ^ 0, and P2^0 there holds 

ui{Q) = Ci, U2{Q) = C2, «3(0) = Ci, W4(0) = -C2, 

^^l(0) = Pl, V2{0) = p2, V3{0) = -pi, V4{0) = P2, 

and 

ui(so) = 0, U2(so) = 0, usiso) ^ 0, U4(so) ¥= 0, 
^^i(so) 7^ 0, V2{so) 0, 1^3(50) = 0, V4{so) = 0, 

then 7(s) extends to a periodic orbit with period Ssq and a symmetry group isomor- 
phic to D4 such that 

wi(3so) 7^ 0, U2(3so) 7^ 0, W3(3so) = 0, U4(3so) = 0, 
wi(3so) = 0, V2{3so) = 0, t;3(3so) 7^ 0, ^4(850) ^ 0, 



and 



and 



ui(5so) = 0, U2(5so) = 0, U3(5so) 7^ 0, U4(5so) 7^ 0, 

^^i(5so) 7^ 0, V2{5so) 7^ 0, V3{5so) = 0, V4{5so) = 0, 

wi(7so) 7^ 0, U2(7so) 7^ 0, U3(7so) = 0, U4(7so) = 0, 

viiTso) = 0, U2(7so) = 0, t;3(7so) 7^ 0, V4{7so) 7^ 0. 
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5. Analytic Existence in the Equal Mass Case 



When m = 1, there is an additional symmetry in the positions of the four masses 
that reduces the PPS4BP with equal masses to the fully symmetric planar four- 
body equal mass problem. We exploit this reduction to prove the existence of a 
symmetric periodic simultaneous binary collision orbit in the equal mass case. 

The additional symmetry is the Ansatz, = xi, x^ = X2, with \x2\ < xi. From 
this it follows that LJ4 = uji, W3 = ijJ2, xi — X2 > 0, X1+X2 > 0. From the canonical 
transformations (1) and (2), we have 



2ul 



= Xi - X2 = 



2ui 



2ui 



= Xi+X2 = 



2ul 



l + ^/2 "'^ -l-|-\/2' l-|-\/2 ' ""^ -l + ^/2' 

Since 2uiU2 = 92= X2 — x\<Q and 2usU4 = 34 = + 0:2 > 0, it follows that 

(7) U2 = -{V2 - l)ui, U4 = {V2- 1)U3. 

Prom the second canonical transformation (2), we have 

Vi = \/2(wi - C02)UI, V2 = -(2 - V2){cOi - i02)ui, 

V3 = \/2(wi + W2)U3, V4 = {2- V2){uJi + UJ2)us. 

These imply that 

(8) V2 = -{V2 - 1)VI, V4 = {V2-1)V3. 

Substitution into the Hamiltonian system of equations with Hamiltonian T (and 
with m = 1) gives 



2^2 2 



U3 = 



4-2^2 



-V3U1, 



(y2-l)(4-2V2) 2 

U2=- '-^ ViU^, 

, (72-l)(4-2x/2) 2 

Ui = V3.U^, 



E' = 0, and 

(4 - 2V2)uivl 



v[=- 



+ 4ui + 



4uiu§ 



5„,2 



^2 



', = -{V2-l) 



(4 - 2V2)uivl 



+ 4ui + 



4miu§ 



,3/2 



+ 2(4 - 2V2)Euiul, 



4uf u| 



,3/2 



+ 2(4 - 2V2)Euiul 



w, = - 



(4 - 2V2)u3vf 



+ 4U3 + 



V<+^3 (Uf + Ui) 



4ufu3 

3/2 



< = (V2 - 1) 



(4 - 2y/2)u3vl 



+ 4U3 + 



4wfu3 



+ 2{4-2V2)Eulu3, 



4ufu3 



4x3/2 



+ 2(4 - 2V2)Eulu3 
t' = {4- 2V2f uluj. 
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Because of Equations (7) and (8), the equations in u'2, u'^, v'2, and v'^ duplicate 
those in u[, Ug, v[, and W3. The Ansatz X4 = Xi, Xs = X2 with \x2\ < Xi, therefore 
leads to the reduced system of equations, 

, 4-2^2 o 

Ul = 1 V1U3, 

, 4-2\/2 2 

U3 = 1 ^^3%, 

E' =0, 



V. =- 



^ V + 4wi + , , , + 2 4 - 2V2)Euiul, 

4 A/4 I ,,4 / 4 , 4\3/2 \ / i ^' 
V ^^1 + % (W| + 1*3) 

, (4-2a/2)u3U? ^ 4u?W3 4u?u^ „ a^x a 9 

^3 = - ^ + + - — + 2(4 - 2x/2)i;.f .3, 

i' =(4 - 2\/2)2u?u^. 
Scale the value of E by 

E 

(9) £; = 



4 - 2V2 ' 



and define 

f = + - 2{ul + ul) - -^$L^ - (4 - 2V2rEulul 

° V -1 + ^ 

It is straight-forward to check that the reduced system of equations satisfies 
, dt , dt . ^ ^ ~, or , df 

Thus the system of reduced equations is Hamiltonian. 

We will simplify the Hamiltonian F by a linear symplectic transformation with 
a multiplier ^ 1. Define new coordinates {Q\,Q2,E,P\,P2,t) by 

(10) Wl = TTTTZ. ^'1 



2V^' ^ 2vVf^' 

2 ' 2J^/2-~l 



(12) E= t = 23/4(^2- 1)3/2 ^. 

^ ' (4-2\/2)2' ^ ^ 

This is a linear symplectic change of coordinates with multiplier 

(13) /z ^ 



25/4^^2 _ 1 

Under this linear symplectic transformation and the accompanying scaling a = s/fi 
of the independent variable s, the Hamiltonian T becomes 



r = 1^ {p!QI + PiQl) - V2{Ql + Ql) - ^f^j^ - EQlQ 
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The reduced system of equations is the Hamiltonian system of equations with 
Hamiltonian T, 

dQi or dPi dV dE dT 

-d^=m^ -da=-m' ' = da=^^^'- 
The function F is a regularized Hamiltonian for the fully symmetric planar four- 
body equal mass problem with the bodies located at {xi,X2), {x2,xi), (— Xi, —X2), 
and {—X2,—xi) (see [15]). On the level set F = 0, the solutions have energy E. 
One regularized simultaneous binary collision occurs when Qi = and Q2 7^ 0, 
for which F = implies = 16\/2, and for which the transformation between 
a-nd xi,X2 implies xi ~ X2 = and Xi + X2 ^ 0. The other regularized 
simultaneous binary collision occurs when Qi 7^ and Q2 = 0, for which F = 
implies P| = 16\/2, and for which the transformation between Q\,Q2 and x\,X2 
implies xi — X2 7^ and xi + X2 = 0. Total collapse occurs when Qi = and 
Q2 = 0, and is the only singularity in F that is not regularized. A solution (5i((t), 
Q2{<j), Pi{o'), P2i(^), 0" € [0,cro], for (To > 0, of the Hamiltonian system equations 
with Hamiltonian F is nonsingular if it avoids total collapse, i.e., + Q2 7^ for 
all cr e [0,ao]. 

The following result is from [15]. The proof of it is a consequence of four equal 
mass bodies starting at (xi,X2), {x2,xi), {—xi,—X2), and {—X2,~xi) with xi = 1 
and X2 = 0, and with the momenta (0, d), {d, 0), (0, —'&), and (— i?, 0) for any i? > 0, 
always having a simultaneous binary collision on the line X2 = xi at a time to > 
continuously depending on «?, such that the cluster velocity x\{tQ) + X2{to) is a 
continuous function of 7?. 

Lemma 5.1. There exists "d > 0, Uq > Q, and a nonsingular solution Qi{a), (32(c), 
Pi(cr), P2{<^), cr e [0, c7o], of the Hamiltonian system of equations with Hamiltonian 
F on the level set F = such that 

(15) Oi(o) = i, g2(o) = i, Pi(o) = -t?, P2(o) = i?, 

(16) E= ^'~'^^^~^ <0, T{a) = £ Ql{y)QUy) dy, 
and 

(17) gi(co) = 0, Q2(co)>0, Pi(ao) - -4(21/4), P2(^(^)=o. 

This Lemma gives the existence of a solution of a boundary value problem for 
the Hamiltonian system of equations with Hamiltonian F. It is this solution whose 
symmetric extension gives a symmetric periodic SBC orbit in the PPS4BP with 

equal masses. 

Theorem 5.2. Fix m = 1. For each E <0, there exists a time-reversible periodic 
regularized SBC orbit 

7(s) = {ui{s),U2{s),U3(s),U4{s),Vi{s),V2{s),V3{s),V4{s)) 

with period T > 0, angular momentum ^4 = 0, and a symmetry group isomorphic 
to D4, for the Hamiltonian system of equations with Hamiltonian F on the level set 
f = such that distinct regularized SBCs occur at s = T/8, 3T/8, 5T/8, 7T/8. This 
periodic orbit corresponds to a symmetric periodic singular orbit 

{Xi{t), X2{t), X3{t), X4{t), UJi{t), U}2{t), W3(i), u;4{t)) 
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with energy E for the PPS4BP with equal masses m = 1 where for all t, 

X4{t) = Xi{t), X-i{t) = X2{t), \x2{t)\ < Xi{t), UJ4{t) = Uli{t), UJsit) = W2 (i) , 

with initial conditions 

xi{0) > 0, a;2(0) = 0, wi(0) = 0, ^2(0) > 0, 

and period 

.T/2 

R= / {u'i{s)+ul{s)){ul{s)+ul{s)) ds, 
Jo 

where for t € [0,R], the only singularities are two distinct SBC's occurring at 
t = R/A and t = 3i?/4. 

Proof. By Lemma 5.1, let Qi{a-), Q2(o"), Pi{o'), P2{o'), € [0,(Jo], be the nonsin- 
gular solution of the Hamiltonian system of equations (14) with Hamiltonian T on 
the level set T = 0, whose properties are given in (15), (16), and (17). Using the 
scaling a — s/fi, set ,so = /vdo. By the linear symplcctic transformation given in 
(10), (11), and (12) with multiplier fj, (as given in (13)), we have 

Q,{s/n) Q2{s/fi) P^js/^i) P2{s/^i) 

Ms) = ' "3(5) = 01 M ' Ms) = — , vais) 



21/4 ' 2V4 ' --^^2V7^' ^''~2VV^1 

By Equations (7) and (8), we have 

-(^/2-l)Ql(s^) (^^-l)Q2(a^) 

^2(3) = ^ , Ui(S) = ^ , 



-VV2-iPi{s/fj,) VV2-lP2{s/^i) 
Ms) = ^ , ^^4(5) = ^ . 

Set 7(s) = {ui{s),U2{s),U3{s),U4{s),vi{s),V2{s),V3{s),Vi{s)), s e [0,so]. Prom the 
Equations (9), (12), (16) in E, E, E, and 1?, we obtain 

E- 3- — <0. 

With this value of E, it follows that the value of F at 7(0) is 0. Set 

Ci = ^>0' C2 = < 0, 



-%') i9 W2 - 1 

Pi = , ^ < 0, P2 = 7, > 0. 



With (5i((t), Q2(c), Pl{<^)^ ^2(0"), cr G [0,(To]; being a nonsingular solution of the 
Hamiltonian system of equations with Hamiltonian F, we have Qi{cf) + Qtif^) 7^ 
for all a G [O,0-o]. Prom this it follows for all s G [0, Sq] that 

Ms = ul{s) - ul{s) + ulis) - ul{s) = {2-V2) [Ql{s/fx) + Ql{s/^i)] ^ 0, 

and 

Ms = 2ui{s)u2{s) - 2u3{s)u4s) = -(2 - V2)[Ql{s/n) + Qlis/fi)] ^ 0. 

These imply that M| + M| ^ and M| + M| ^ for all s G [0,so]. Thus the 
function 7(5) is a nonsingular solution of the Hamiltonian system of equations (3) 
with Hamiltonian F on the level set F = that satisfies 

Ul(0)=Cl, «2(0) = C2 W3(0)=Cl, W4(0) = -C2, 
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Vi{0) = pi, V2{0) = P2, V3{0)=-pi, V4{0) = P2, 
Ul(so)=0, U2(so)=0, Us{so)>0, U4{so)>0, 
Vi{so)<0, V2{so)>0, W3(so) = 0, t;4(so)=0. 

By Lemma 4.1, the solution 7(3) extends to a T = 8so periodic solution, call it 
7(5), with a D4 symmetry group generated by the symmetries Sp and Sg, and four 
distinct regularized simultaneous binary collisions at s = sq, 3so, 5so, 7so, for which 

(so) + ul{so) = 0, ul{so) + ul{so) ^ 0, 
uj{3so) + ul{3so) ^ 0, uliSso) + ul{3so) = 0, 

(5so) + ul{5so) = 0, ul{5so) + uI{5sq) 7^ 0, 

u?(7so) + ul{7so) ^ 0, ul{7so) + uI{7sq) = 0. 

Since Qfia) + Q|(cr) ^ for all a e [0,ao], it follows that + Mi(s))(Mi(s) + 

^K'S)) 7^ foi' s e [0,T] except s = so,3so,5so,7sq. The regularizing change of 
time (4), 

£ = (ulis) + ulis)){ulis) + ul{s)), 

defines t as an invertible different iable function of s, i.e., t = 6{s) with ^(0) = 

and 9'{s) = when s = {2k + 1)sq for fc G Z. The symmetry Sp satisfies 
Sf7{s) = 7(s + 2so) and —7(5) = Sp^{s) = 7(s + 4so)- The symmetry Sg satisfies 
Sg7{s) = 7(2so — s), and so 7(5) has a time-reversing symmetry. The angular 
momentum (5) for 7(5) at s = is 

A = ^ [ - V1U2 + V2U1 - V3U4, + V4U3\ = P2C1 - P1C2 = 0. 

The extended solution 7(s) gives a singular symmetric solution 

Z{t) = {xi{t),X2it),X3{t),X4{t),UJi{t),UJ2{t),u;3{t),U)4{t)), 

of the PPS4BP with m = 1. Under the Ansatz, the components of z{t) satisfies 
X4{t) = xi{t), xs{t) = X2{t), \x2{t)\ <xi{t), uji{t) = uJi{t), uj3{t) = ui2{t), where 

^ uljs) - uljs) + uljs) - ul{s) 

X2{t) = Ui{s)u2{s) + ■U3(s)u4(s), 

_ Vl{s)Ui{s) - V2{s)u2{s) V3{s)u3{s) - V4{s)U4{s) 

_ Vl{s)u2{s) + ^2(5)^1(5) V3is)u4{s) + ^4(5)^3(5) 

for s = d~^{t). The components of the extended solution 7(5) satisfy 113(0) = ui(0), 

W4(0) = -U2{0), «l(0)U2(0) < 0, |W2(0)| = (72- l)|Ui(0)|, V3{0) = -Vi{0), V4{0) = 
V2{0), 



(0)1x2(0) + 1^2(0)1x1(0) = P1C2 + P2C1 = ^^^4^ > 0. 



21/4 



and 



W2(0)ii2(0) - i;i(0)wi(0) = P2C2 - PiCi 



1 I? 



21/4 

From Lemma 2.1, it follows that a;i(0) > 0, a;2(0) = 0, wi(0) = 0, and ^2(0) > 0. Set 
R = 6{T/2). Since 7(4so) = -7(0), it follows that xi{R) = xi{0), X2{R) = X2{<S), 
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wi(i?) = u!i{0), and C02{R) = ^2(0). Thus the singular symmetric solution z{t) has 
period R. By the construction of the extension of 7(5) given in Lemma 4.1, there 

holds 

(fe+l)so /■So 

{ul{s) + ul{s)){ul{s) + ui{s)) ds= (u?(s) + ul{s)){ul{s) + ul{s)) ds 

kso J 

for all fc = 1, . . . , 7. This implies that R/A = e{{k + l)so) - 0{kso) for all k = 

0, 1, . . . , 7. The first regularized SBC for 7(5) occurs at s = s^, and this corresponds 
to t = 9{so) = R/4:. The next regularized SBC for 7(5) occurs at s = 3so, and this 
corresponds to 

t = e{3so) = m^o) - 0i2so)) + (0(2so) - 0{so)) + 0{so) = | + | + | = ^. 

Similarly, the regularized SBCs for 7(5) occurring at s = 5so,7so correspond to 
t = 5i?/4, 7i?/4. Hence, for t G [0,R], the periodic solution z{t) has SBCs as its 
only singularities, and these occur at t = -R/4, 3JR/4. 

For a fixed but arbitrary e > 0, the value of €~'^E is a fixed but arbitrary negative 
real number. By Lemma 3.1, the scaled extended solution 

7e(s) = (eui(es), eU2{es), eus{es), eU4{es),Vi{es),V2{es), V3{es), v^ies)) 

is a periodic solution for the Hamiltonian system of equations (3) with Hamiltonian 
r on the level set F = 0, having period e~^T and energy e~^E <0. By an argument 
similar to above, 7c (s) satisfies the required conditions. □ 

6. Numerical Estimates in the Equal Mass Case 

In the equal mass case, there is by Theorem 5.2 a time-reversible periodic orbit 
7(5) for the Hamiltonian system of equations with Hamiltonian F on the level set 
F = with period T. The components ?i2(s), u^i-s), U4{s), vi{s), V2{s), 

^3(5), V4{s) of 7(5) satisfy Equations (7) and (8). The D4 symmetry group of 
7(5) is generated by 5f7(s) = 7(s + T/4) and Sc'yis) = 7(^/4 — s). Under the 
linear symplectic transformation (10), (11), and (12) with multiplier fi (as given 
by (13)), this gives a periodic orbit Qi{a), Q2{cr), Pi{cr), P2{'y) of the Hamiltonian 
system of equations (14) with Hamiltonian F on the level set F = 0, which by 
Lemma 5.1 satisfies Qi(0) = 1, 02(0) = 1, Pi(0) = -t?, P2(0) = t?, (3i(cro) = 0, 
<32(o-o) > 0, Pi(c7o) = -4(21/4)^ and ^2(0-0) = for some ctq > and > 0. In [3], 
we numerically estimated 

£70 = 1.62047369909693, 1? = 2.57486992651942. 

The period of this periodic orbit is 8cro ~ 12.96378959 and its energy \s E k, 
—2.999682732. From the linear symplectic transformation (10), (11), and (12), with 
multiplier ji and from Equations (7) and (8), we have for the values of components 
of 7(0) the exact 

wi(0) = W3(0) = 2-V4, -U2{Q) = Ui{Q) = {^^2- l)2-^'\ 

and the estimates 

vi{Q) = -U3(0) = , ss -2.000382939, 

2V\/2- 1 



t;2(0) = ^^4(0) = "^^"^ ^ w 0.8285857433. 
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Figure 2. The ui, U2, U3, U4, vi, V2, fs, V4 coordinates of the 
symmetric periodic SBC orbit 7e(s) for m — 1 and an e > 1. 

Since a = s/n, the period of 7(5) is T « 8.469003682. From Equation (9), the 
value of the energy for 7(5) is E -5.120778733. 

The components of the scaled periodic orbit 7e(s) for an e > 1, shown in Figure 
2, satisfy the symmetries Spjeis) = 7e(s + T'e/4) and Sc'^eis) ~ — s). We 

choose e so that = 27r, and check the linear stability of 7e(s) (see Theorem 
3.3). Using a Runge-Kutta order 4 algorithm with a fixed time step of 27r/50000, 
we computed X^{2'k). Two of the eigenvalues of X^{2-k) are 1 by Theorem 3.2. 
Numerical estimates of the remaining eigenvalues of X^iiir) are 

-0.9888731375 ± 0.1487612779i, 
-0.9973584383 ± 0.07263708002i, 
0.9999060579 ± 0.01370676220i, 

that have modulus one. Thus numerically, the periodic orbits 7e(s) are linearly 
stable for all e > 0. The first complex conjugate pair of eigenvalues for X^{2'k) 
matches the complex conjugate pair of characteristic multipliers for the periodic 
orbit (3i((t), Q2{<j)i -Pi(o'), ^2(0') of the Hamiltonian system of equations (14) with 
Hamiltonian F, corresponding to 71 (s) = 7(s), where we computed [3] the real 
part of the complex conjugate pair of modulus one to be —0.9888840619. Because, 
by a lengthly computation, JV^r(7(:(s)) is not block diagonal, the last two com- 
plex conjugate pairs of eigenvalues of X^(2'k) are not repeats of the characteristic 
multipliers of the periodic orbit Qi((t), Q2(o'), A(o'), ^2(o')- 

Figure 1 illustrates the curves in the physical plane that the four equal masses 
follow in the linearly stable SBC orbit z{t) — {xi{t), X2{t), X3{t), X4{t), LLi{t), uj2{t), 
uJait), uJ4{t)) of the PPS4BP, corresponding to 7<:(s) with e = 1/(2 - v^). The 
initial conditions for z{t) are 

a;i(0) =0:4(0) = 1, X2(0) =X3(0) = 0, 

cJi(O) = CJ4(0) = 0, a;2(0) = ^3(0) « 1.287434964. 
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The value of e = l/v2 — \71 here for the scahng is determined by the equation 
xi(0) = e'^{ui{0) — ""2(0)) coming from Lemma 2.1 and the canonical transforma- 
tions (1) and (2) applied to the scaled periodic solution 7e(s), together with the 
initial condition a;i(0) = 1, where t = corresponds to s = 0. The value of the 
Hamiltonian H along z{t) is e'^E w -2.999682732. 

7. Numerical Estimates in the Unequal Mass Case 

We numerically continue to < m < 1 the time-reversible periodic regularized 
simultaneous binary collision orbit for the Hamiltonian system of equations with 
Hamiltonian T on the level set F = 0. The analytic existence of this orbit is given by 
Theorem 5.2 for m = 1. We also numerically estimate the monodromy matrices and 
their eigenvalues for these continued time-reversible periodic regularized simultane- 
ous binary collision orbits for < m < 1. For m = 1, we assume by Lemma 3.1 that 
the time-reversible periodic regularized simultaneous binary collision orbit 7(5; 1) 
with the D4 symmetry has period T = 27r and energy E w —2.818584789. Recall 
that the 1)4 symmetry group of j{s; 1) is generated by SF'y{s; 1) = 7(,s + tt/2; 1) 
and <S'(37(s; 1) = 7(7r/2 — s; 1). We assume without loss of generality by Lemma 
3.1 and Lemma 4.1, that the continued time- reversible periodic regularized simul- 
taneous binary collision orbit 7(5; m) for < m < 1, has period T ~ 27r, energy 
E{m), and a D4 symmetry group generated by 5^7(5;™) = 7(5 + 7r/2;TO) and 
Sgj{s; m) = 7(7r/2 — s; m). We assume E{m) depends continuously on < to < 1, 
with E{1) = E. We shift the regularized time variable s to s -|- 7r/4, so that for 

7(s; m) =(ui(s; m), U2{s\ m),U3{s; m), ^4(5; to), 
vi{s; m),V2{s; m),V3{s; to), V4{s; to)), 

the value s = now corresponds to the first simultaneous binary collision, i.e., 
■ui(0;to) = 0, M2(0;to) = 0, U3(0;to) ^ 0, U4(0;to) 5^ 0, t;i(0;TO) ^ 0, V2{0;m) ^ 0, 
113(0; to) = 0, and V4{0,m) = 0. We approximate the continued time-reversible 
periodic orbits 7(5; to), < to < 1, by the trigonometric polynomials, 

n n 

ui{s; to) = Qi sin((2z — l)s), U2{s; to) = — bi sin((2i — l)s), 

i=l i=l 

U3{s;m) = Ui{s — 7r/2; to), U4(s;to) = U2{s + 7r/2;TO), 

n n 

vi{s;in) = cos((2i — l)s), V2{s;m) = — ^(iiCos((2i — l)s), 

i=l i=l 

W3(,s; to) = vi {s — 7r/2; to), ^4(5; to) = V2{s + 7r/2; to), 

for a positive integer n and constants Qi, bi, Ci, di, i = 1, . . . ,n, that are assumed to 
be continuous functions of to. The presence of odd positive integer frequencies 2i — l 
in these trigonometric polynomials is to ensure that the period functions defined 
by them have the D4 symmetry group generated by Sp and Sg- So in particular, 
the periodic orbits in the continuation are time-reversible for all < to < 1. An 
numerical estimate of the periodic solution 7(s; to) is found through the variational 
approach of minimizing the functional 

/■27r 

L= ||7'(s; to) - JVf(7(s; to))|| ds 
Jo 
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Figure 3. The estimated values of 1/3(0; m), U4{0;m), vi{0;m), 
and ^2(0; m) for the nonzero components of 7(0; m) for < m < 1. 



over the space M of coefficients a^, 6^, d^, z = 1, . . . , n, for an appropriate choice 
of n. The use of trigonometric polynomials for numerically approximating periodic 
solutions is a classic approach (see, for example, Simo [25]). 

The numerical algorithm for finding a trigonometric polynomial that approxi- 
mates the periodic solution 7(5; m), < m < 1, proceeds in two steps. First, we 
consider a guess for the set of values for Ui, bi, Ci, di, i = 1, . . . , n, as well as a guess 
for E{m). Starting with a reasonably low number of terms, n = b, we let a numer- 
ical minimization algorithm (in this case, MATLAB's fminunc) find the minimum 
of L near the starting guess. Then we add an additional non-zero term to each 
of the trigonometric polynomials, and the minimizing solution from the previous 
iteration is used as a starting guess for the next iteration. This process continues 
until we reach n = 10. Second, since the Hamiltonian equations with Hamiltonian 
r requires a specified value of Eim) to compute L, we need to make certain that 
we are getting a good estimate of E{m). We evaluate F at the point 7(7r/4;m), 
away from simultaneous binary collisions. If this value is not sufficiently close to 
(within about 5 x 1G~^° of 0), we adjust E{m) using the bisection method in a 
small interval about the initial guess of E{m) until r(7(7r/4; m)) is sufficiently close 
to 0, re-minimizing the trigonometric polynomial approximation of 7(5; m) for each 
new choice of E{m). 

This numerical method only works well if we have a good initial guesses for a^, 
Ci, di, I = 1, . . . , 5, and E{m) for some value of m. This we have when m = 1. 
We use our estimate of 7(5; 1) and E{1) to provide the initial guesses for Oi, bi, 
Ci, di, i = 1, . . . , 5, and E{m) for m = 0.99. The numerical algorithm produces a 
trigonometric polynomial approximation of 7(5; 0.99) and an estimate of E{0.99), 
which then provide the initial guesses for a^, bi, Ci, di, i = 1, . . . ,5, and E(m) for 
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Figure 4. The energy E{m) of the periodic orbits 7(5; m) for 
< TO < 1. 




Figure 5. The value of f(7(7r/4;TO)) for the optimized trigono- 
metric polynomial approximation of 7(3; to) over < m < 1. 

TO = 0.98. We continue decreasing m by 0.01 and using the numerical algorithm 
until we reach m = 0.01. In Figure 3, we plot the graphs of the numerical estimates 
of U3(0;to), M4(0;m), t;i(0;m), and U2(0;to). Notice that ui(0;to) and W2(0;to) 
satisfy Equation (6), as is expected from the regularization of the simultaneous 
binary collisions. In Figure 4, we plot the graph of the numerical estimate of E{m). 
In Figure 5, we graph of value of r(7(7r/4; m)) over O < m < 1. 

We use the optimized trigonometric approximations of 7(5; to), < to < 1, to 
estimate 7(— 7r/4;TO), which are the conditions corresponding to s = before the 
shift of s. For each < m < 1, the conditions 7(— 7r/4; to) satisfy Lemma 2.1, and 
so they corresponds to the initial conditions of interest as given in the Introduction. 

We numerically estimate the monodromy matrices X{s; to) and their eigenvalues 
for 7(5; to) with < TO < 1. We used a Runge-Kutta order 4 algorithm coded in 
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Figure 6. The maximum modulus of the eight eigenvalues of 
X{s;m) for 0.2 < m < 1. 

MATLAB with a fixed time step of 27r/50000 at increments of 0.01 in m. In Figure 6 
we graph the maximum modulus of the eight eigenvalues of X{s; m) over 0.2 < m < 
1. The maximum modulus of the eigenvalues of X{s]m) for < m < 0.2 becomes 
large as m — >■ 0, and are therefore not included in Figure 6. The maximum modulus 
for the eight eigenvalues shows that 7(5; m) is linearly unstable for < m < 0.53. 
The error in the numerical estimation of the repeated eigenvalue 1 of X{3]m) is 
large enough to account for the graph in Figure 6 not lying at 1 over 0.54 < m < 1. 
However, for 0.54 < m < 1, the other six eigenvalues come in three distinct complex 
conjugate pairs that are each of modulus one and located away from ±1. This shows 
that 7(5; m) is linearly stable for 0.54 < m < 1. 

Acknowledgements. We thank the referees for their valuable comments that 
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